QUANTUM SYSTEMS OBEYING A 
GENERALIZED EXCLUSION-INCLUSION 

PRINCIPLE 

A.M. Scarfone 1 

Dipartimento di Fisica - Politecnico di Torino 
Corso Duca degli Abruzzi 24, 10129 Torino, Italy; 
Istituto Nazionale di Fisica della Materia - Unita del Politecnico di Torino, 
Corso Duca degli Abruzzi 24, 10129 Torino, Italy. 

29, february 2000 



Abstract 

The aim of this work is to describe, at the quantum level, a many body system 
obeying to a generalized Exclusion-Inclusion Principle (EIP) originated by collec- 
tive effects, the dynamics, in mean field approximation, being ruled by a nonlinear 
Schrodinger equation. The EIP takes its origin from a nonlinear kinetic equation 
where the nonlinearities describe interactions of different physical nature. 
The method starts from the study of the kinetic behavior of many particle system 
which results to be nonlinear because of the interaction among the particles and 
introduces an effective nonlinear potential £/eip which permits us to simulate the 
true interactions governing the dynamics of the system. The power of the method is 
tested in the case of spatially homogeneous classical iV-particles system. Its kinetic 
in the momenta space is described by a Markofnan process. A judicious general- 
ization of the particle current, permits us to obtain at the equilibrium, a statistical 
distribution interpolating in a continuous way the well known quantum statistics 
(Fermi-Dirac or Bose-Einstein). 

Systems with statistical behavior interpolating between the Bose-Einstein and the 
Fermi-Dirac were introduced fifty years ago. Up to now many have been the at- 
tempts to generalize quantum statistics. As we will discuss in the first chapter 
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many of these generalized statistics can be obtained by means of EIP by an ap- 
propriate generalization of the many particle current expression. The result of this 
approach is a nonlinear Schrodinger equation with a complex nonlinearity. Only the 
imaginary part of this nonlinearity is fixed by EIP, while the real one, is obtained by 
the requirement that the system is canonical. This permits us to include together 
within EIP other interactions of different physical origin. 

Systems obeying to EIP can be used to describe a wide range of physical situations. 
In condensed matter we can find many topics where EIP concepts can be applied 
like, for instance, in superconductivity to describe the formation of vortex-like exci- 
tations or in superfluidity to describe the formation of Bose-Einstein condensates. 

The plane of the work is the following. In the first chapter, after an introduction, 
we explain the origin of EIP. In the same chapter we summarize the fundamental 
mathematical tools used in the thesis: variational principle, Lagrangian and Hamil- 
tonian formalism for systems with infinity degrees of freedom, symmetries and con- 
servation quantities. In chapter II, EIP is obtained from an appropriate deformation 
of the quantum current. From all the possible deformations of the current we chose 
the simplest one to study the mean properties of systems obeying to EIP. By means 
of variational principle we obtain a canonical nonlinear Schrodinger equation. We 
study the Lagrangian and Hamiltonian formulation both in the standard ip repre- 
sentation than in the hydrodynamic one. In chapter III we look at the symmetries 
of the nonlinear Schrodinger equation obeying to EIP. By means of Nother theorem 
we obtain and discuss conserved quantities associated with the symmetries of EIP 
potential and show that the systems obeying to EIP are conservative. Only internal 
forces are introduced in the system by the presence of EIP potential. In the same 
chapter we consider an important class of nonlinear gauge transformations. These 
transformations are obtained in order to linearize the continuity equation. As a con- 
sequence, the transformed evolution equation will be again a nonlinear Schrodinger 
equation containing a purely real nonlinearity. In chapter IV EIP-Schr6dinger equa- 
tion is coupled in a minimal way to an abelian gauge field. The dynamics of the 
gauge fields is described by the Maxwell-Chern-Simons Lagrangian. In presence of 
Chern-Simons coupling the evolution of the system is developed in (2+1) space-time 
dimension so that the model can be used to describe planar phenomena. This situ- 
ation, for instance, is realized in the layer between two semiconductors or between 
superconductor and normal conductor and in particular in High-T c superconduc- 
tors where Chern-Simons coupling is used to describe the anyonic behavior of the 
excitations. We show that in this case, the anyonic statistic behavior ascribed to 
the system by the presence of Chern-Simons Lagrangian is not destroyed by the 
presence of EIP potential. The following two chapters are devoted to the study 
of special solutions of evolution equations of the system. In chapter V we study 
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the solitary waves in neutral systems. Here we obtain the differential equation for 
the shape of the solitons, and study its mean properties. Using the nonlinear gauge 
transformations obtained in chapter III we deduce the effective potential responsible 
of the formation of the solitons. Applications on the Bose-Einstein condensation are 
given. In chapter VI we consider charge, static, self-dual vortex-like solutions where 
the dynamics of gauge fields is ruled by the Chern-Simons Lagrangian. The shapes 
of the vortices as well as its electric and magnetic fields are obtained by numerical 
integration of appropriate differential equations. The physical properties of this so- 
lution are matched to the corresponding solutions known in literature where EIP is 
absent. Conclusions are reported in the last chapter. 

PACS: 03.65.-w, 05.20.Dd, 11.15.-q 
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Chapter I 



Introduction to EIP 



1.1 Considerations on quantum 
many body systems 

Macroscopic systems are characterized by a high number of degrees of freedom that 
makes impossible to determine exactly their evolution. This is due to many reasons 
like, for instance, the large number of evolution equations that must be solved (one 
for each degree of freedom), the incomplete knowledge of the initial conditions, the 
approximate knowledge of the physical interaction and so on. It also true that for 
complex systems it is more important to have information about physical observable 
that represent mean quantities characterizing the system as a whole. 
We deal with the evolution equations describing the dynamics of the mean values 
of these observable. Because of the high number of degrees of freedom, it is more 
convenient to study complex systems in the context of field theories. Fields carry 
infinite degrees of freedom and their dynamics is ruled out by partial differential 
equations. The self-interactions introduced by the collective effects of the many 
particle systems are in general described by nonlinear terms in the evolution equa- 
tions. 

An important question in the topics of many body system is the statistic behavior of 
its constituents. This question is even more important for a quantum system where 
concept of statistics plays an important role. 

In this thesis we deal with a nonlinear theory describing in the mean field approxi- 
mation an interacting many body system where the nonlinear interactions are intro- 
duced starting from the kinetics of the system. Therefore the topics of the nonlinear 
partial differential equations and the statistic behavior of a quantum many body 
system are important subjects of this thesis. Let us clarify these two arguments. 
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We begin with a review about quantum statistics. Since the early days of quan- 
tum mechanics, it was clear that, from the principle of undistinguishability, the 
statistical behavior of a collection of identical particles must be different from the 
classical one. The particle statistics determines the structure of the many body 
wave functions, that turn out to be completely symmetric under permutations of 
identical particles named bosons which obey to the Bose-Einstein statistics or com- 
pletely antisymmetric under permutation of identical particles, now called fermions 
and obeying to the Fermi-Dirac statistics. 

The statistics obeyed by fermions and bosons have many important implications in 
quantum mechanics. For instance, the Pauli exclusion principle, which gives null 
probability for two fermions to be in the same quantum state, takes its origin from 
the antisymmetric fermionic wave functions; it has consequences like the quality of 
emitted spectrum of the atoms or the stability of compact objects likes white dwarfs 
or neutron stars. On the contrary, the totally symmetric bosonic wave function has 
the effect of enhancing the probability for the bosons to occupy a quantum state if 
this one is already occupied. This effect is responsible of phenomena like supercon- 
ductivity or superfluidity. 

In the early fifty Green [1, 2] found that the principles of quantum mechanics allow 
two kinds of statistics called para-Bose and para-Fermi statistics [3]. The parastatis- 
tics of order p are defined as the identical particles statistics in three dimensional 
space under the restriction of a possible number of particles in the symmetric or 
antisymmetric state for the para-Fermi and para-Bose respectively. The case with 
order p — 1 corresponds to the ordinary Fermi and Bose statistics. The different 
cases can be described by trilinear commutation relations among the creation and 
annihilation operators [4]. The parastatistics was applied to subnuclear components 
like quarks [5] to solve, for example, the puzzle in the quantum number of the bari- 
onic resonance A ++ . 

The case of a particle described by a value of p not an integer has been studied in 
Ref. [6, 7], in order to take into account small violations of Pauli exclusion principle 
or Bose statistics. In Ref [8] this particle was called paronic. However, the corre- 
sponding quantum field theories for such a particle turns out to have negative norm 
states and, as a consequence, are not acceptable [9, 10]. This saga culminates with 
a recent study of infinite statistics without assumptions on the parameter p [11, 12], 
all representations of the symmetric group can occur. The particles obeying this 
type of statistics are called quons [11, 13, 14]. The quonic statistics are described 
by the g-deformed bilinear commutation relations: 

aa) — qa' a = 1 , (1.1.1) 

where q is a C-number with \q\ < 1 and according to the Fredenhagen theorem [10], 
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cannot be embedded in the local algebra of observable. Differently from the paronic 
statistics the quonic one has positive definite squared norms for state vectors but 
notwithstanding Greenberg succeeded only in the nonrelativistic quantum theory 
due to the locality problem for the infinity statistics. 

Differently from the quonic statististics, the concept of g-deformed oscillators or 
q-oscillators [15, 16, 17] takes its origin from the concept of quantum groups [18]. 
g-oscillator, in spite of classic oscillators obeys to the g-deformed unitary algebra 
SU q (N): 

aa^Tqa j a = q TN , (1.1.2) 

with q G C and N is the number operator functions of the annihilation and creation 
operator a and a* . In (1.1.2) the sign — or + is reported to the g-bosons or q- 
fermions, respectively. The g-oscillators can formally be defined in any dimensional 
space but then violate the fundamental axioms of quantum field theories in terms 
of the relation between spin and statistics. 

The situation changes radically in many body systems, the dynamics being con- 
fined in two spatial dimensions. As it was discussed by Wilczek [19, 20], the statistic 
behavior of a system depends on the property of interchange of identical particles 
and is related to the topological property of the configuration space of a collection of 
identical particles. In more than two dimensions only two possibilities are present. 
Here the fundamental group is that of permutation which has two one-dimensional 
representations corresponding to completely symmetric (Bosons) or antisymmetric 
(Fermions) wave functions. Differently, in two spatial dimensions the permutation 
group is replaced by the braid group [21]: the spin is not quantized in integer 
or half integer value and particles obey any statistics interpolating between the 
Bose and Fermi ones. These particles are called anyons. Successively in Ref. [22] 
were achieved similar conclusions using a completely different method based on the 
study of the unitary representation of the current algebra and diffeomorfism group. 
Anyons can occur in many physical applications in those condensed matter systems 
that can be effectively regarded as two dimensional. For example, anyons can occur 
in fractional quantum Hall effect where collective excitations have been identified as 
localized quasi-particles of fractional charges, fractional spin and fractional statistics 
[23, 24, 25] and also they can occur in high temperature superconductors recently 
discovered [26]. 

Is not too difficult to show in a heuristic way how anyonic statistics born in two 
spatial dimensions. Let 2) be the function describing two identical particles 
and assume that when we move particle 2 around particle 1 by an angle A ip the 
wave function changes as: 

VK1, 2) ^ 2) = e- A ^(l, 2) , (1.1.3) 
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acquiring a phase factor depending on a statistical parameter v. Now we look at 
the exchange of the two particles. This can be realized in two ways (see figure 1.1): 
moves particle 2 around particle 1 by an angle Atp = n or at the opposite side, by 
an angle A(p = —tt. In the two cases the wave function acquires an extra phase 
exp(±7i v). It is easy to recognize that in more than two spatial dimensions the two 
paths are topologically equivalent, so identifying the two transformations we have 
the relation: 

e inu = e -i*»^ (1.1.4) 

which is true only if v — 0, 1 modulo 2, corresponding, respectively, to the well 
known bosonic and fermionc statistics. In more than two dimensions there are no 
other possibilities. Differently, in two dimensions we can not deform with continu- 
ity the paths, one into the other. They are topologically and physically distinct 
operations. The Eq. (1.1.4) does not necessarily hold any more and the statistical 
parameter v can be not shrunken to take the value or 1. The anyonic statistics 
was obtained in a more rigorous way by Y. -S. Wu [27]. For a review on anyons see 
for example Ref. [28]. 

We have emphasized the anyonic statistics because in chapter IV and VI we discuss 
a particle system obeying an exclusion-inclusion principle, where the matter field is 
coupled to an abelian gauge field whose dynamics is described by means of Chern- 
Simons Lagrangian. As it was discussed in Ref. [19, 29, 30] the presence of the 
Chern-Simons gauge field confers the anyonic behavior to the system. 

Another definition of generalized quantum statistics has been formulated by Hal- 
dane [31, 32], is based on the rate of the number of the available states in a system 
of fixed size decreasing as more particles are added to it. This statistics is called 
exclusion statistics. The statistics of Haldane is formulated without any reference to 
spatial dimensions of the system. In his formulation of exclusion statistics, Haldane 
defines a generalized Pauli exclusion principle introducing the dimension of the 
Hilbert space for single particle states as a finite and extensive quantity that depends 
on the number N of particles contained in the system. The exclusion principle im- 
plies that the number of available single particle states decreases as the occupational 
number increases 

Ad N = -gAN , (1.1.5) 

where g is the parameter that characterizes the complete or partial action of the 
exclusion principle and makes possible the interpolation between the Bose-Einstein 
(g = 0) and the Fermi-Dirac (g = 1) statistics. The relevance of the Haldane 
statistics is in its implication in the fractional quantum Hall effect and anyonic 
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physics, in the Calogero-Sutherland model [33, 34, 35] and in the Luttinger model 
[36, 37]. 

It is known [38] that the effects due to the statistics are imposed by the Pauli 
exclusion principle to a system of free fermions and can be simulated by a repulsive 
potential in the coordinate space. Analogously, free bosons can be submitted to 
an attractive potential. We refer to it as a statistical potential. The statistical 
potential will be a nonlinear function of the fields describing the system and its 
spatial derivatives. 

Several nonlinear Schrodinger equations (NLSEs) have been studied in the past 
and recently, they are commonly used in many different fields of research in physics. 
The cubic equation [39], for instance, with the nonlinearity proportional to ±|?/>| 2 , 
has been used to study the dynamical evolution of a boson gas with ^-function 
pairwise repulsion or attraction, responsible of its anyonic-like behavior [40]. Re- 
cently, this equation has been used to describe the Bose-Einstein condensation 
[41, 42, 43, 44, 45] and the dynamics of two-dimensional radiating vortices [46]. 
The nonlinear term \ip\ 2 appears also in the Ginzburg-Landau model of the super- 
conductivity [47], a phenomenon investigated also by means of the Eckhaus equation 
which is a NLSE with a nonlinearity of the type |^| 2 +a \ifj\ 4 [48]. The same equation 
appears in superfluidity, where the properties of a gas of bosons interacting via a 
two-body attractive and three-body repulsive S function inter-particle potential are 
investigated [49, 50]. The Eckhaus equation can describe nonlinear waves in optical 
fibers with a " normal" dependence of the refractive index on the light intensity [46] . 
Another important example where nonlinearities in the Schrodinger equation induce 
a statistical behavior is given by Schrodinger-Chern-Simons theory. In Ref. [51] it 
was shown that the gauge fields can be expressed as functions of the matter fields 
and therefore can be eliminated from the initial equation by transforming it into a 
highly nonlinear Schrodinger equation which describe the same anyonic system. 
In literature we can find NLSEs with complex and derivative type nonlinearities 
involving the quantities (Vp) 2 , Ap, j ■ Vp, V • j [52, 53, 54] as, for instance, in the 
Doebner-Goldin equation associated with a certain unitary group representation and 
describing irreversible and dissipative quantum systems. NLSEs with nonlinearities 
involving the quantity j have been also introduced to study planar systems of par- 
ticles with anyon statistics [55]. 

We will see in chapter II that the potential introduced by the EIP is complex and 
derivative in the field and ifi*. 

Solution spectrum of nonlinear partial differential equations is much rich than the 
linear one. Generally, the solutions of a nonlinear PDE can be split in two distinct 
classes: in the first we have the solutions that can be obtained with the perturbative 
method. This is possible if the coupling constant of the nonlinearity is small. In the 
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other class we find the nonperturbative solutions, obtained integrating directly the 
nonlinear PDE. This solutions depend on the coupling constant of the nonlinearity 
and are divergent in the zero limit. So, it is not possible to go continuously from the 
solutions of perturbative class to the nonperturbative class. In this one we find the 
soliton solutions. They are wave packets that propagate freely asymptotically not 
changing their shape and velocity also after a collision. These were discovered in the 
last century by Russel, but their interest in physics was emphasized only thirty years 
ago. Solitons are solutions of nonlinear partial differential equation. For instance, if 
we take into account the dispersion relation of the linear Schrodinger equation we 
have: 

h 2 k 2 

nuJ = — , (1.1.6) 
and looking for the group velocity v g = dk/dui: 

< 117 » 

which depends on the number wave k. This means that each component in the 
wave packet propagates with a different velocity and therefore the modulation of 
the shape of the packet is not maintained. Responsible of this dispersion is the 
Laplacian operator appearing in the kinetic term. To avoid this effect a confining 
potential is required which makes the evolution equation nonlinear. 
Soliton solutions were found in many fields, for example light impulse in the wave 
guide [56, 57, 58, 59], propagation in electric circuits and plasma waves [60]. Their 
particle behavior makes them interesting in the physics of elementary particles [61] 
like for instance the t'Hooft-Polyakov monopole [62, 63]. 

The question of the research of the solutions of nonlinear differential equations is one 
of the most important topics of mathematical physics of the last years. One of the 
most powerful is the inverse scattering method [64] . It is a canonical transformation 
in the action-angle fields in which the Hamiltonian of the system appears to be diag- 
onalized. This transformation is a Fourier transformation plus a Laplace one. When 
an evolution equation is solved with this method, the Hamiltonian spectrum results 
decomposed in a continuous part given by the perturbative solutions plus a dis- 
crete contribution given by the soliton solutions. Equations solved with the inverse 
scattering method are named ^-integrable because their solutions are expressed as 
function on a spectral parameter. The more easy S-integrable equation is of course 
the Schrodinger equation which is solved by means of Fourier transformation. 
Finally a new class of integrable equations have made their appearance recently. 
These are called C-integrable because can be linearized by means of change of de- 
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pendent variables. In this category we can find for example the Burger equation [65] 
or the Ekhaus equation [66]. 

1.2 What is EIP 

We present the generalized Exclusion-Inclusion Principle (EIP) in the configuration 
space. For a rigorous introduction of the EIP we remind to the reference [67, 68, 69, 
70, 71, 72]. In the first part of this section the discussion is keep at a classic level. 
The EIP takes its origin from a classical nonlinear kinetic that takes into account 
inhibition or enhancement of the particle transition probabilities in the phase space. 
We start by considering a Markofnan process in a D-dimensional phase space. If 
we identify the state of the system by a vector in the phase space and setting 
7r(t, u — > v) the transition probability from the state u to the state v, the evolution 
equation for the distribution function n(t, v) can be written as 

dU ^ ^ = / [n(t, «->«)- n(t, v -> «)] d D u . (1.2.1) 

The exclusion-inclusion principle is introduced into the classical transition probabil- 
ities by means of an inhibition or an enhancement factor. 

In Ref. [68] it was postulated the following expression of the transition probability: 

7r(t, v — > u) = r(t, v, v — u) (f)[n(t, v)] ip[n(t, u)] , (1-2-2) 

where r(t, v, v — u) is the transition rate, <f>[n(t, v)] is a function depending on the 
occupational distribution at the initial state v and ip[n(t, u)\ depends on the arrival 
state. The function <p{n) must obey the condition 0(0) = because the transition 
probability is equal to zero if the initial state is empty. Furthermore, the function 
ip(n) must obey the condition ip(0) = 1 because, if the arrival state is empty, the 
transition probability is not modified. The classical linear case is obtained if we 
chose 0(n) = n and ip(n) = 1. 

In the following, for reason of simplicity, we consider an one-dimensional space, in 
the first neighbor interaction approximation. 

If we consider an infinitesimal transition from the state v to the state v + dv, the 
transition rate can be defined as [68] 

r(t, v, ±dv) dv 2 = D(t, v) ± )-J(t, v) dv , (1.2.3) 

where J(t, v) and D(t, v) are the drift and diffusion coefficients, respectively. In the 
simple case in which in a time interval dt only one transition is allowed, Eq. (1.2.1) 
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becomes 



dn(t, v) 

at 



n(t, v — dv — > v) + ir(t, v + dv — > v) 
ir(t, v — > v — dv) — ir(t, v — > v + dv) . 



(1.2.4) 



The (1.2.4) is a balance equation between the particle coming in v ± dv — > v in the 
site v with respect to that coming out v — > v ± dv. It describes in the continuum 
limit the discrete Markoffian process reported in figure 1.2. 

Expanding the r.h.s. of Eq. (1.2.4) in powers of dv, up to second order, in the limit 
dv — > and taking into account the transition rate defined in Eq. (1.2.3), we obtain 
the following generalized, nonlinear Fokker-Planck equation: 



dn(t, v) 



d_ 

dv 



J(t, v) + 



dD(t, v) 
dv 



+ D(t,v) U( n )^|^-0(n) 



(f)(n) ip{n) 
dijj(n)\ 



dv rv ' dv 

This is a continuity equation for the distribution function n = n(t, v) 

dn(t, v) dj(t, v, n) 



1.2.5) 



dt dv 
where the particle current j = j(t, v, n) is given by: 



0, 



(1.2.6) 



J 



' . d D(t, v)\ . . . . ^. , / . ,d(b{n) , . .dib(n) 
J (t, v) + <t>{n) i){n) + D(t, v) *p(n)-^ - 0(n)- 



dv 



dv 



dv 



1(1.2.7) 



When the system reaches the equilibrium configuration, the current (1.2.7) must be 
equal to zero. In particular, if we consider Brownian particles, the drift and the 
diffusion coefficients are given by: 



J 



7 v 



D = 



7 



(5 m 



(1.2.8) 



with 7 a dimensional constant and (3 = 1/ksT where ks is the Boltzmann constant. 
In this situation, when j = Eq. (1.2.7) can be viewed as a first order differential 
equation with solution: 



ip(n) 



= e 



[1.2.9) 
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where e = (3 (E — fi) and E = mv 2 /2 is the kinetic energy. The integration constant 
/i is the chemical potential which can be evaluated by fixing the number of particles 
of the system. 

We show now how we can select the functions <f>(n) and ip(n) in order to obtain 
the equilibrium distribution of some of the statistical distributions described in the 
previous section. 
If we make the choice: 

4>(n) = n , ip{n) = 1 + ku , (1.2.10) 

from (1.2.9) we obtain the following stationary distribution: 

n=—^—. (1.2.11) 

Here we recognize the well know Maxwell-Boltzmann (MB), Bose-Einstein (BE) 
and Fermi-Dirac (FD) distributions when we fix the value of the statistic parameter 
k — 0, ±1 respectively. Moreover for — 1 < k < 1 we have fractional distributions 
interpolating between the BE and FD ones. 
As a second example we pose: 

(j>(n) — n (1 — gn)^~ [1 + (1 — g) n] * , 



(1.2.12) 

ij(n) = (l-gn) 1 ¥ll + (l-g) n} 1 ^ , (1.2.13) 

and obtain the distribution of Haldane's statistics: 

ne e = (1 -gn) 9 [l + (1 -g)^ 1 ' 9 , (1.2.14) 

obtained in Ref. [32]. The statistical parameter g is defined in (1.1.5). Finally if we 
pose: 

<P(n) = [n] q , = [1 + Hn (1.2.15) 
where [x\ q is defined as: 

[x] q = £^-, (1.2.16) 
we obtain the g-oscillator distribution: 

e-= Sinh[r ^ + " n)] , (1.2.17) 
sinh(r)n) 
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with r] = log q. The parameter a — ±1 permits us to describe bosonic or fermionic 
g-oscillator respectively. 

We can give other examples but it is now clear how to obtain different statistical 
distributions starting from a nonlinear kinetic. 

In this thesis we shall study, in mean field approximation, a canonical quantum 
system obeying to an exclusion-inclusion principle obtained by choosing <f>(n) = n 
and ip{n) = 1 + ku. 

This choice is made because it is the most easy to treat and, as we have shown be- 
fore, permits us to simulate a system with an equilibrium distribution interpolating 
between the BE and FD distribution. 



Now we introduce EIP in a quantum system. 
Let us start by looking at an exclusion-inclusion principle in the configuration space. 
We consider the classical stochastic Marcoffian process in a 3-dimensional space 
described by the following forward nonlinear Fokker-Planck equation (FPE): 

^ + V.j«=0, (1.2.18) 
at 

with 

=u (+) p(l + K P ) -DVp , (1.2.19) 

where p = p(t, x) is the occupational number or particle distribution in the config- 
uration space, u( + ) is the forward velocity and k E IR. 

Beside the forward FPE (1.2.19) we must also consider the backward one with cur- 
rent: 

j(-)= u (-)p(l + «p) + DVp. (1.2.20) 
The semisum of the forward and backward FPEs will give us 

-^ + V[vp(l + K p)] = 0, (1.2.21) 



where v is the current velocity given by: 

1 



v — - 
2 



M W + «(-) 



(1.2.22) 



We stress one more on the meaning of the current. We define the transition proba- 
bility from the site x to x' as n(t, x — > x') = r(t, x, x') p(t, x) [1 + k p(t, x')\ with 
r(t, x, x') the transition rate and the choice <p{p) = p, ip = 1 + k p. The transition 
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probability depends on the particle population p(t, x) of the starting point x and 
also on the population p(t, x') of the arrival point x' . If k > the n(t, x — > x') 
introduces an inclusion principle. In fact the population at the arrival point x' stim- 
ulates the transition and the transition probability increases linearly with p(t, x'). 
In the case k < the 7i(t, x — > x') takes into account the Pauli exclusion principle. 
If the arrival point x' is empty p(t, x') = 0, the 7r(£, x — > £c') depends only on the 
population of the starting point. If the arrival site is populated < p{t, x') < p max 
the transition is inhibited. The range of values the parameter k can assume is 
bounded by the condition that 7r(t, x — > x') be real and positive as the r(t, x, x'). 
Then we conclude that k is limited from below by the condition k > —l/p max . 

In the Nelson picture [73, 74, 75], a quantum system can be viewed as a stochastic 
process where the particles are subjected to a Brownian diffusion with a coefficient 
D = h/2m. The quantum system is described by the complex wave function if) = 
ip(t, x) and the quantity p = \ip\ 2 is interpreted as the particle probability density. 
We make the ansatz: 

i> = p 1 ' 2 exp (is) , (1.2.23) 

where the phase S is related to the velocity v as [75, 76]: 

V S 

v = — . (1.2.24) 

m 

Using Eqs. (1.2.23) and (1.2.24) it is immediate to obtain from Eq. (1.2.21) the 
expression for the quantum current: 

j = -i^- (i + Kp ) (r vip-ip vr) , (1.2.25) 

which we can rewrite also as: 

vs 

3 = p(l + K P ) , (1.2.26) 

m 

which is our starting assumption. 

In Ref. [77], by using Nelson stochastic quantization method, it was obtained the 
following NLSE: 

^ |^ = A i, + A(p, j)1> + V1>, (1-2.27) 

where: 

A(p,j) = W(p)+iW(p,j) , (1.2.28) 
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is the complex nonlinearity introduced by the EIP, and the expression of the imag- 
inary part is: 



j) 




) 



(1.2.29) 



The real part W(p) depends drastically on the quantization method. In the picture 
of stochastic quantization it was obtained the form: 



It can be shown that the system described by Eqs. (1.2.27), (1.2.29) and (1.2.30) 
admits a continuity equation (1.2.18) with the current given by Eq. (1.2.25). 
We define in the following the quantum system obeying to EIP as a system whose 
dynamics is described by a NLSE admitting a continuity equation with the quantum 
current gave by Eq. (1.2.25). 

We conclude this section by considering some examples where EIP can be use- 
fully applied. In nuclear physics the correlation effects between pairs of nucleons, 
viewed as fermions, are quite relevant in the interpretation of experimental results. 
Similarly, the interactions among bosons are relevant in various nuclear modes (su- 
perfluid model, interacting boson model, mean field boson approximation) and allow 
the explanation of many collective nuclear properties. The interaction among the 
fermionic valence nucleons outside the core produces pairs of correlated nucleons that 
can be approximated as particles with a behavior intermediate between fermionic 
and bosonic ones. This nuclear state (quasideuteron state) can be viewed as a par- 
ticles system that obeys to EIP. Recently, it was studied a semiclassical model of 
photofission in the quasideuteron energy region [78] . We described the quasideuteron 
state as a mixture of fermion and boson states, with a good agreement of our cal- 
culated photofission cross sections of several heavy nuclei and experimental results. 
Another example is the Bose-Einstein condensation. The condensation originates 
from an attraction of statistical nature (Bose-Einstein statistics) among the parti- 
cles. In several papers the Bose-Einstein condensation is studied by means of a cubic 
NLSE [41, 42] which describes in mean field approximation an attractive interaction 
between two bodies. In place of the cubic and simplest interaction, other interac- 
tions can be considered as, for instance, the one introduced by EIP to simulate an 
attraction among the particles. 

In condensed matter we can consider the problem of the hopping transport on a 
lattice of ionic conductors: two ions having the same charge cannot occupy the 
same site due to their natural electrostatic repulsion; also the motion of the couple 



W(p) = K 



ft 2 Ap 2 - k p 



+ 



(V P ) 2 




4m |_1 + np 2p (1 + up) 2 



12 



electron- hole in a semiconductor can be described by means of EIP. In fact, while 
electrons and hole are fermions, together can be considered excited states behaving 
differently from that of a fermion or a boson. 



1.3 Mathematical background 

In this section we summarize the canonical method used in the study of systems 
with infinite degrees of freedom whose dynamics is described by a partial differential 
equations. 

Let us start with some definitions and notations. Let M to be a complex smooth 
manifold with dimension D mapped by the coordinates Xj with % — 1, • • • , D. Let 
T be the algebra of the functionals on M — > M of the type: T 3 P = JV d D x. 
We call the quantity V functional density. Let t/j(t, x) and ip*(t, x) be two fields on 
M 2 , we set \1/ = (■?/>, ip*) the complex 2- vector field on M x M. We assume that the 
2-vector field \l/ vanishes quickly on the boundary of M. 

We consider now a non relativistic canonical quantum system described by the 
Lagrangian density £[ip] = C(d t ip*, d t ip, [ip*], [ip]) € T depending on the scalar field 
ip E M and its derivatives with d t = d/dt. Here and in the following we adopt the 
notation in square bracket to indicate the dependence from the spatial derivative of 
any order. Moreover, for our purpose, we deal with the case in which the spatial 
derivative in the Lagrangian density are introduced by V operator. Thus, the 
following notations means: 

C([a]) = C(a, Va, V 2 a, V 3 a, • • •) . (1.3.1) 

Posing the action functional: 

A = J Cd D xdt , (1.3.2) 

the evolution equations for the fields ip and ip* can be obtained by a variational 
principle [79]: 

5A = 0, (1.3.3) 
where the variation of a functional F e T is a 2-vector defined as: 

(1.3.4) 




2 Because the theory developed in this thesis is nonrelativistic, here and after t play the role of 
a parameter 
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The functional derivative can be defined by means of Euler operator: 



which it is given by: 



E, = EE(-^) J ^, (1-3-6) 



n=0 J ^ J 



where the second sum is extended over multi-indices J = (j t , j±, • • • , jo) with < 
ji < n, % — t, 1, • • • , D and jt + J2 ji — n - Eq. (1.3.3) are the Euler-Lagrange 
equations for the field ip* and ip respectively. From Eq. (1.3.3) we obtain: 

£ ( - 1),V "8(^-^ = - (1-3.7) 
and its conjugate equation. 

It is possible to show that the Euler operator satisfies the following property: 

E(H + VC) = 0, (1-3-8) 

with B, C G T . Therefore, Lagrangian density which are total derivatives does not 
gives contribute to the evolution equations. It are named null Lagrangian. 
In this thesis, we will consider canonical systems described by the following class of 
Lagrangian density in (3+1) dimensions: 

c = 4 { r Tt " ^ TFt) " w|2 " um ' m - y ^^' (1-3-9) 

where V describes an external potential and £/([■?/>*], [ip]) is a real nonlinear potential. 
Using the Lagrangian (1.3.9) we obtain the following NLSE for the field ip: 

which is a derivative NLSE. 

In a canonical theory, we can express the motion equations also in the Hamilton 
formalism. This requires that functional space M is a Poisson manifold, i.e. a 
.D-dimensional smooth manifold equipped with the Poisson brackets. This one is 
defined as: 

{F, G} = J (5FV5G) d D x , (1.3.11) 
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where F, G £ T . V is the Hamiltonian linear operator acting on T which might 
depends by ip and its derivative. 

The Poisson brackets must satisfy the following properties: 

• Linear 

• Skew-symmetric 

• Jacobi identity 

This property is reflected on the definition of V. The first property is evident from 
the linearity of T> and the definition of Poisson brackets, the second require T> to be 
skew-adjoint: = —V. Finally, the third condition is: 

{ {P, Q}, R} + cyclic permutations = . (1.3.12) 

This last relation is true if V is a constant quantity not dependent on the fields ip 
and ip*. 

Eqs. (1.3.7) can be expressed in the Hamiltonian formalism if a functional T 3 H = 
J Hd D x called Hamiltonian function exists (H is the Hamiltonia density) so that: 

— = V5H, (1.3.13) 
dt y ' 

which is equivalent to the following expression in terms of Poisson brackets: 



dt 



H} . (1.3.14) 



In general, given a functional F e T describing a physical observable of a system 
with Hamiltonian H, its time evolution can be written as: 

d F 8 F 

£ = + (1.3,5, 

We remark en passant that for a system with finite dimension, the Poisson brack- 
ets must satisfy one more condition: it must be a derivative operator which means 
to satisfy the Leibnitz rule {A B, C} = A{B, C} + {A, C} B. For infinitely di- 
mensional systems this condition is not imposed because the multiplication between 
elements of T is not well defined. In fact, it gives two functionals A, B E J 7 , their 
product is not expressible as integral of a density functional i.e. the relation: 

J Ad D x J Bd D x = J Cd D x , for same density C (1.3.16) 
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is not generally satisfied. 

Let us now introduce the fields ir^ and n m , canonically conjugate momenta of the 
fields ip and ip*, respectively and define: 



dC 



dC 



7T, 



7T, 



d(d t r) 



[1.3.17) 



We can write the Hamiltonian density Ti.([ip], [71^ ]) related to the Lagrangian density 
by the Legendre transformation [80]): 



n 



dip 



dip* 



dt + V dt 



(1.3.18) 



and the evolution of the system is described by the equations: 

SH 



d^ 
~dt 

dt 



5 n , 



1.3.19) 



5H_ 



It is easy to see that Eqs. (1.3.19) are equivalent to Eq. (1.3.13) if we make the 
choice for the Hamiltonian operator: 



V 



1 
-1 



1.3.20) 



which is skew-hermit ian and satisfies the Jacobi identity because it has constant 
entry. 

With this definition for the operator V the expression (1.3.15) becomes: 



dt r 

u J 



5 ip(z) 5 ir^ (z) 5 ip(z) 5 tt^ (z) 
8F(x) 8H(y) 5H(y) 5F(x 



d D z 



6iI)*(z)6tt (z) 5iP*(z)5tt (z) 



d D z + 



dF 
~dt 



(1.3.21) 



which can be simplified following the Dirac procedure [81] (see sections 2.1 and 4.2). 
An important fact, in canonical theories, follows from the Nother theorem [82], which 
states the link between symmetries for the Lagrangian (or Hamiltonian) and the 
conserved physical quantities related to it. We resume briefly the general method. 
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The Nother theorem says that for each one-parameter symmetry group of the system, 
there is a physical observable related to it: 

Q = J J°d D x , (1.3.22) 

which is conserved as a consequence of the continuity equation: 

8 7° d V 

-w+ikr*' (L3 - 23) 

where i — 1, • • • , D (sum on i is assumed) and J % are the components of a D- vector 
describing the flux density associated to the J° . 

In fact, let Sip be the variation produced by the action of a group transformation 
on ip, it generates a symmetry for the system if the action (1.3.2) does not change. 
This implies for the Lagrangian: 

5C = d v f v , (1.3.24) 

which says that it might change for a total divergence (null Lagrangian). To give 
an example, we consider the easy case in which only the first derivative of the fields 
are present in the Lagrangian density. Computing the variation in C we obtain: 

and after integration by part, taking into account the motion equations (1.2.21) we 
arrive at Eq. (1.3.23) where: 

We remember that the J v is not univocally defined. In fact, we can add to the 
expression (1.3.26) the gradient of a skew-symmetric 2-tensor: 

SF" — > lf v — \J V + <9 M , (1.3.27) 

with 

TV » = _ T v» _ (1.3.28) 

As a consequence of the skew-symmetry the new vector J v satisfies the same con- 
tinuity equation (1.3.23) and the conserved quantities are left unchanged by the 
substitution (1.3.27) if the field goes to zero on the boundary of M. 
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Chapter II 



Canonical Systems Obeying to the 
Exclusion-Inclusion Principle 



In this chapter we introduce a class of canonical nonlinear Schrodinger equations 
obeying to a generalized inclusion-exclusion principle (EIP). This is accomplished 
trough an opportune deformation of the expression of the quantum current of the 
matter field ip. The Lagrangian and Hamiltonian structure are studied both in the 
^-representation and in the hydrodynamic one. The approach used to include the 
EIP in the evolution equation does not determine in a unique way the form of the 
nonlinear potential. In fact, we show that in the nonlinear potential U EIP , which is 
a complex functional, its real part is not determinable from kinetic considerations. 
We impose the canonicity of the system in order to determine the real part of the 
nonlinear potential. The final result is different from those obtained by performing 
other quantization method. Moreover, the canonicity required leave us the possi- 
bility to include an arbitrary real nonlinear potential U[p\. This arbitrarity can be 
used to introduce other interactions acting on the system simultaneously with the 
EIP. Finally we study the mean property of the system. In particular we analyze 
the Eherenfest relations for observables like energy, linear and angular momentum. 

2.1 Canonical systems 

We start from the general expression of a nonlinear Schrodinger Lagrangian in (D + 
1)— dimension of the form: 




dip* 



) 



%■ 



2 



£ = i- 



dt 



2 m 



wr- w*], m-vr^P, (2.1.1) 
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where V = (d/dxi, d/dx 2 , ...d/dxo) is the gradient operator in Z)-dimension, 
H is a constant with the dimension of an action that we identify with the Planck 
constant and m is the mass parameter. 

The first two terms in the Lagrangian density are the same encountered in the stan- 
dard linear quantum description, the quantity V is a potential describing external 
interaction. U([ip*], [if)]) is the nonlinear term which we assume to be an analytic 
smooth functional of the fields ip, tp* and their spatial derivatives. We assume U 
real to make the system not dissipative. 
Using the Lagrangian we introduce the action: 



A 



j Cd D xdt . (2.1.2) 



Applying the Euler operator E^„, defined in Eq. (1.3.6), to Eq. (2.1.2), we obtain 
the following NLSE for the field ip: 

">^ = -&^ + W m ' %m + Vi, < (2 ' L3) 

where (A = d 2 /dx 2 + d 2 /dx% + • • • + d 2 /dxp) is the D-dimensional Laplacian 
operator. 

To obtain the expression of the nonlinear potential U([ip*], [ip]) in order to include 
the EIP in the system it is convenient to consider the Bohm-Madelung representation 
for the wave function ijj: 



ip(t, 1 x) = pit, x) 1/2 exp 



j:S(t, x) 



(2.1.4) 



The hydrodynamic fields, density of particles p and phase S [83] are related with ip 
by means of: 

P=\M 2 , (2-1.5) 

The nonlinear potential U([p], [S]) becomes now a real functional of the fields p, S 
and their spatial derivatives. 

By taking into account the Leibnitz rule for the functional derivative: 

5 5 p 5 S S 5 5 . h , S , n _ „. 



5ip* 6ip*6p 6tp*5S op 2p T SS 
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Eq. (2.1.3) is rewritten in the form: 



.8^ h 2 A , 



j- p U([ P ], [S]) 



i ■ h 

2p 



Introducing the nonlinear quantities: 

W([ P ], [S]) 



T-U{\p\ y [S]), 
bp 

MM, W) = j- p js u{[pl [s]h 



ip + Vip .(2.1.8) 

(2.1.9) 
(2.1.10) 



the NLSE (2.1.8) becomes: 



ih^ = -^Arl, + W{\p], [S])ij + iW([p}, [S])4> + V^. (2.1.11) 

Using this expression it is easy to obtain the continuity equation of the system. By 
taking the product of (2.1.11) times ip* and subtract the complex conjugate form, 
we arrive to the expression: 



where 



dp 
dt 



Jo 



h 



i h 



2m 



(2.1.12) 



(2.1.13) 



is the standard quantum current density. 

Following the discussion made in section 1.2, we require that the system described by 
(2.1.11) admits the following expression for the quantum current [cfr. Eq. (1.2.25)]: 



j = (l + K P )j , 
and that the continuity equation (2.1.12) becomes: 

d p 



dt 



+ V • j = . 



(2.1.14) 



(2.1.15) 



This can be accomplished if, taking into account expression (2.1.14) and (2.1.12), 
we select for the term W the form: 



2p 



2p \l + K Pj 



(2.1.16) 
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which has the same form of Eq. (1.2.29) obtained in Ref. [77] using the stochastic 
quantization method. Eq. (2.1.16) can be expressed with the hydrodynamic fields 
as: 



W([rt,[S]) = -«Av(^) , (2.1.17) 

and taking into account the dependence of W from the nonlinear potential U([p], [S]) 
given by Eq. (2.1.10) we obtain: 

U{\p\AS]) = K { ^p* + U{[ P \) , (2.1.18) 

Zi lit 

which is defined modulo a real functional of the field p and its spatial derivative. 
We note that starting from the continuity equation only the dependence of U([p], [S] ) 
from the phase can be determinate whilst the dependence from the field p is not 
fixed which means that an arbitrary real quantity [S]) functional of p, can 

be added compatibly with the EIP, without removing the canonicity of the system. 
Hereinafter we call " EIP potential" the quantity C/eip(p, S) given by: 

U ElP (p, S) = K^^p 2 . (2.1.19) 
2m 

This potential depends on p and on V5. Inserting Eq. (2.1.18) in Eq. (2.1.9) we 
obtain for the real part of the nonlinearity the expression: 

W(p, S) = k- (j^—) 2 + F([p}) , (2.1.20) 

p \l+Kp) 

with 

F(\P\) = 4-U(\p\) ■ (2-1.21) 
bp 

We have obtained the following result: 

Starting from the expression of the quantum current j appearing in the continuity 
equation, it is possible to deduce a NLSE compatible with it. This NLSE generally 
contains a complex nonlinearity. Only its imaginary part is fixed while the real one 
required an additional constraints. If we require that the system obeying the EIP is 
canonical, we obtain for the real part the quantity W(p, S) gives by Eq. (2.1.20) 
that is defined modulo a real functional F([p\). 

In Ref. [77] the iter of the stochastic quantization leads to a NLSE with the same 
expression of W but a different form for W [see Eq. (1.2.30)]. 
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Returning again to the fundamental fields ip and ■0*, the expression of the po- 
tential (2.1.19) takes the form: 

u E1P m, m) = -K-?-(rvi>-i>vr) 2 , (2.1.22) 

so that we can write the most general form of the Lagrangian of a system obeying 
the EIP as: 

2 \ at at J 2m 8m 

- U(U\ 2 ])-Vr^, (2.1.23) 

where t/QI^I 2 ]) is a real arbitrary smooth functional which might depend on p and 
its spatial derivatives. 

In this thesis, if not otherwise specified, the arbitrary potential U(p) is an analytic 
nonderivative functional of p only. This potential can be used to describe other 
interactions in the system. By an appropriately choice of its form, the Lagrangian 
(2.1.1) can be used to describe different physical systems in presence of collective 
interactions and obeying to the EIP (see chapter V). 
Using Eq.(2.1.3) we obtain the following NLSE: 

f) l fi 2 fi 2 

ih^= - — a^-k — M>* v^j -i> v^j*) w 

ot 2m 2m 

h 2 

- K — v{ip*vi)-i)Vi)*)i) + F( P )i) + vi), (2.1.24) 

where F(p) is now given by F(p) = dU(p)/d p. 

Eq. (2.1.24) can also be written using the expression of the current (2.1.14) and the 
particle density p in the form: 

iH It = ~^- A ^ + A ^> 3) i> + F{p) ip + Vip , (2.1.25) 
ot 2m 

where the complex nonlinearity A(p, j) is given by: 

Mp,J)=n^( T J-)'-in^-v( T ^-) . (2.1.26) 

p \l + KpJ 2p \1 + KpJ 

The quantum system described by the Lagrangian density (2.1.1) is canonical. 
This can be verified defining the fields ir^ and n^,, canonically conjugated to the 
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field ip and ip*, by means of the relations (1.3.17): 

■ h i* 



7T* 

4> 



■ k / 
—l — W. 

2 y 



(2.1.27) 
(2.1.28) 



It is well known that tt^ and ir , are proportional to the fields ip* and ip, so that, while 
in the Lagrangian formalism ip and ip* are independent fields, in the Hamiltonian 
formalism they are canonically conjugated. Following Ref.[81], Eqs. (2.1.27) and 
(2.1.28) give rise to the primary constrains: 



(2.1.29) 
(2.1.30) 



Performing the Legendre transformation (1.3.18), it is easy to see that the Hamil- 
tonian density can be written as: 



n 



2m 



h 2 



' m 



(V>* - ^ W) 2 + £/(V>* i)) + Vi)*i) . (2.1.31) 



Let us introduce now the Poisson brackets between two functionals: 



{/(*), s(y)} 



<$/(*) Sg(y) 5 f(y) 6 g(x) 



+ 



I 



6 ij>(z) 6 n^z) 5n^(z)5ip(z 
8f(x) Sg(y) 5 f(y) 5g(x 



d D z 



_5r{z)5^{z) 5^{z)5r{z)\ d z - (2 - L32) 

The second class primary derivative (2.1.29) and (2.1.30) satisfy the relation: 

{6(t, x), &t, y)} = -th5^ D \x - y) , (2.1.33) 
and can be accommodated by the introduction of the Dirac brackets: 

{f(x),g(y)} D = {f(x),g(y)}+ l -l{f(x),^(z)}{Uz),g(y)}d D z 

" k I {f{xh Uz)} {Uzh 9{y)} d ° z ■ (2 - L34) 

Expression (2.1.32) can be simplified if one solves the derivative (2.1.29), (2.1.30) 
for tt^ and tt^ and treats / and g as functionals of tp and ip* only. A straightforward 
calculation gives: 



{/(*), 9(v)}d = 



I 



Sf(x) Sg(y) 5g(y) 5 f(x) 



5ip(z) 5ip*(z) 6i/j(z)6il)*(z) 



d D z .(2.1.35) 
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Using this expression we can obtain the evolution equations for the fields ip and ip* 
in the Poisson formalism: 



Ch\)_ 

~dt 
~~dt 

which give the Hamiltonian equations 



= {^,H}, (2.1.36) 
= {r,H}, (2.1.37) 



that are the Schrodinger equations for the fields ip and ip*, respectively. It can 
be easily verified that the chose equations are equal to Eq. (1.3.13) introduced in 
section 1.3 with the Hamiltonian operator given by (1.3.20). 



2.2 Hydro dynamic formulation 

In this section we describe the NLSE with EIP in the hydrodynamic representation. 
Its utility will be seen in chapter V, where we study explicit solutions of the model. 
In particular we make use of the hydrodynamic formulation in order to obtain the 
solitary wave solutions of the system. 

It is well known [83] that a quantum system can be seen as a Madelung-like 
fluid [84] described by means of the fields p and S trough a coupled system of 
differential equations. The velocity field associated to the fluid is related to the 
phase of ip from the relation v = V S/m and its dynamics is described by means of 
a Hamilton- J acobi like equation in presence of a nonlinear term function of p which 
was called quantum potential. Let us remember now that the dynamical equation 
for the phase is only formally similar to the Hamilton- Jacobi equation. In fact, in 
the classical mechanics, the Hamilton- Jacobi equation describes the evolution of the 
Hamilton principal function S that is the generator of a canonical transformation 
and is identified a posteriori with the action of the system: S — J Ldt. On the other 
hand, differently from the Bohm picture, the phase S of the field ip is a dynamical 
field canonically conjugate to the density p and is not a generator of canonical 
t r ansf or mat ions . 

By introducing Eq. (2.1.4) in Eq. (2.1.24) and separating the real part from the 
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imaginary one, we obtain a coupled system in p and S: 
at 2m m 2m y/p 



9 P 

/t +v 



p(l + Kp) 



m 



. (2.2.2) 



Equation (2.2.1) is a Hamilton- J acobi type equation, where the fourth term is the 
quantum potential [83]. The third term is the real part W(p, j) of the term intro- 
duced by the EIP potential (2.1.19), as can be verified taking in mind the expression 
of the current j given by Eq. (2.1.14), that in the new field becomes: 

VS 

3 = p(l + K P ). (2.2.3) 

m 

Finally the last two terms in Eq. (2.2.1) are the extra nonlinearity and the external 
potential. We can see that the imaginary part W(p, j) of the term introduced by 
the EIP potential does not appear in the Hamilton- Jacobi equation (2.2.1). 
It is easy to recognize that Eq. (2.2.2) is the continuity equation in the p-S repre- 
sentation. 

The Hamilton- Jacobi equation (2.2.1) and the continuity equation (2.2.2) can be 
derived from a variational principle, making use of the Euler operator (1.3.6) written 
in the p and S representation, obtained from Eq. (1.3.6) with the substitution ip — > p 
and ip — > S respectively. In fact, starting from the Lagrangian density C which is 
given by: 

„_ as p _(v^ p _ K (vsi p2 _ pvpr_ 0(p) _ Vpi (2 . 2 . 4) 



dt 2m 2m 8m p 

and applied to the action A[p, S] defined by (1.3.2) the Euler operator: 

E p (A) = , E S (A) = , (2.2.5) 

we obtain the Eqs. (2.2.1) and (2.2.2). In the Lagrangian density (2.2.4) the 
first three terms are the same as occurring in the linear Schrodinger equation, 
while the nonlinear contribution is given by the forth and six the term, where 
the fourth is the potential introduced by the EIP. The term -/i 2 (Vp) 2 /(8m/)) 
in Eq. (2.2.4) is responsible of the presence of the quantum potential in Eq. (2.2.1), 
U g = —h 2 Ap 1 / 2 /(2mp 1 / 2 ). Notwithstanding, in literature the quantum potential is 
referred to the nonlinear term that appears in the Hamilton-Jacobi equation. 
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(2.2.8) 
(2.2.9) 



Now we introduce the Hamiltonian procedure. Let us show the results without 
describing the Dirac procedure. 

The momentum n s , canonically conjugate to S, is given by Eq. (1.3.17), that now 
becomes: 

vr s = -p . (2.2.6) 

Moreover, we have 7r p = ik/2. Therefore, ir s is proportional to p, while ir p is a 
constant; the number of degrees of freedom is the same in both the Lagrangian and 
Hamiltonian formalism. 

The Hamiltonian density, function of the canonically conjugate fields S and — p, can 
be deduced taking into account Eq. (1.3.18) and (2.2.4): 

* - ^+^ + ^ + ^- 

The Hamilton- Jacobi and the continuity equations take the form: 

d_S SH 
~dt ~ ~Tp~ ' 
dp _ 5H 

~dt~Ts ' 

The same equations, in the Poisson formalism, can be rewritten as: 

|f = {S, H} , (2.2.10) 

|f = {P, H} . (2.2.11) 

The evolution equations (2.1.38), (2.1.39) or (2.1.36), (2.1.37) deduced from H(ip, ip*) 
have the same form of the equations (2.2.8), (2.2.9) or (2.2.10), (2.2.11) deduced from 
H(p, S). Then, according to a well-established procedure, we can relate the fields 
■ip-ip* to the fields S-p by means of a canonical transformation [80]. The equations 
of motion in the S-p representation will be used in chapter V to study particular 
soliton solutions of Eq. (2.1.25) that preserve their shapes in the time. As we will 
show in chapter V, we are able to decouple the system of equations (2.2.8), (2.2.9) or 
equivalently Eqs. (2.2.10), (2.2.11) obtaining a differential equation in the variable 
p only, whose solutions define the solitons of the systems with EIP. 

In conclusion of this section we study the most simplest solutions of the NLSE 
(2.1.24) in presence of the nonlinearity introduced by EIP only, i.e. when we neglect 
the potentials involving U and V . These are the planar waves: 

ip(t, -x) — A exp [i(k -x — ut)] , (2.2.12) 
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where k is the wave number and A a complex constant. After inserting Eq. (2.2.12) 
in (2.1.24) we obtain the following dispersion relation for the planar waves: 



hk 2 



(1 + 2k\A\ 2 ) , 



(2.2.13) 



U) = 



2m 



which reproduces the well known dispersion relation of the Schrodinger equation 
when the EIP is switched off (k — > 0). Let us remark that in presence of EIP, the 
angular frequency depends also on the amplitude of the wave. A this point it is 
important to remember that, differently from the linear theory, in the NLSE the 
normalization constant is not arbitrary. Generally, the amplitude of a solution is 
related to the other parameter of the system. For instance, the cubic NLSE [60] 
admits soliton solutions whose velocity is related to the amplitude. This amplitude 
increases when the soliton velocity decrease. 

We remember also that because of the nonlinear nature of the systems the more 
general solution can not be expressible as superposition of planar waves how it 
is made in the linear case, by using the Fourier transformation method. More 
sophisticated mathematical tools, like, for instance, the inverse scattering method 
[64] are needed to found the general solutions. In this thesis we do not develop this 
procedure. 

From Eq. (2.1.31), it appears that the system may be unstable in the case k < 
-1/\A\ 2 . In fact, the energy of the modes with wave number k is: 



which must be positive. This imposes the condition for the repulsive systems 



according to the exclusion principle in the configuration space [77] 

2.3 Physical observable 

Now we study the time evolution of the average of the most important physical 
observable that describe the system in presence of EIP. We will identify the motion 
constants of the system. The results of this section will be obtained again in the next 
chapter from the analysis of the symmetries satisfied by the system. The section is 
organized in to parts: in the first we describe the general results. The mathematical 
proofs are collected in the second part. 




2 m 



(2.2.14) 




(2.2.15) 
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2.3.1 Ehrenfest relations 

Let us assume that the nonlinear potential U(p) and the field ib vanish at infinity so 
that the surface terms can be disregarded. Moreover we assume that the potential 
U (p) depends on the space and on the time only through the field p. In this section 
we assume D = 3. 

To obtain the Ehrenfest relations of the system obeying to Eq. (2.1.25) we start, at 
first, with the definition of average of an Hermitian operator A — A': 

<A>= J i\)*Ai\)d 3 x . (2.3.1) 

Hereinafter we normalize the system as: 

N = [ib*ibd 3 x , (2.3.2) 



where N is an integer, so that the field p assumes the meaning of a density of prob- 
ability of position of a iV-body system. Of course the normalization is maintained 
in the time thanks to the continuity equation (2.1.15): 

d JL = . (2.3.3) 
dt 

This means that the number of particles of the system do not change during the 
evolution. Using Eq. (1.3.21) it is easy to obtain the following relationship for the 
time evolution of the average of A: 



d -7 % 
— < A >= - 

dt h 



— Aib-ib* A — 
Sib Sip* 



OA 

d 3 x+ < — > . (2.3.4) 
at 



Here the last term takes eventually into account explicit time dependence of the 
functional A. Let us call O the operator in the r.h.s. of the NLSE (2.1.25) that can 
be rewritten in the form: 

ih^ = Oib ) (2.3.5) 

with 

O = H + W(p, j) + i W(p, j) + F(p) , (2.3.6) 

where H = (—h 2 /2m) A + V(x) is the Hamiltonian operator of the linear theory. 
We can write the relation (2.3.4) in the following form: 

j- < A >= % - < [Red, A] > +± < {ImO, A} > + < ^ > , (2.3.7) 
dt n n at 
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where the symbols [•, •] and {•, •} stand for the commutator and the anticommutator, 
respectively (see also Ref. [53] where an example of the Ehrenfest relations in a 
nonlinear Schrodinger equation with complex potential is discussed). 

After setting A = x c in Eq. (2.3.7), where x c = x/N, we obtain the Ehrenfest 
relationship for the time evolution of the center of mass frame: 



d . h 

— < x c >= —i — 

dt 2mN 



J(1 + k P ) (^*V^-^W) d 3 x . (2.3.8) 



It is worth remarking that the EIP introduces the additional quantity < k, p P > 
which is equal to Kmpj/(\ + k p). Note also that the right hand side of Eq. (2.3.8) 
can be written as: 

d If 

— < x c >= — / j d 3 x , (2.3.9) 
dt N J y ' 

which appear formally the same as in the standard linear quantum mechanics. 
A second relation is obtained setting A = P: 

^- <P>=- j i,* W^d 3 x- j ^j*V F^d 3 x . (2.3.10) 

Equation (2.3.10) can be regarded as the second law of the dynamics [85, 86]. The 
dynamics of the mean value of the momentum is governed by an effective potential 
given by the sum of the external potential V and of the nonlinearity F(p). The 
EIP potential does not affect, on the average, the dynamics of the system because, 
due to their particular form, the terms W and W satisfy the relation < [W, V] — 
% {W, V} >= 0. For the most frequent nonlinearity F(p) generally appearing in the 
nonlinear Schrodinger equations the last term can be dropped and the Newtonian 
behavior is restored [85] (see the next section for the proof of this statement). On 
the contrary, other dynamical equations, like the sine-Gordon equation, seem to 
show a different behavior with respect to the Newtonian one. In this case the reason 
is in the kink-like solution of this equation that do not vanish on the boundary at 
infinity so that the surface terms can not be neglected. 

Than, setting A — L, where L is the angular momentum operator whose com- 
ponents are L i = e^Xj Pk, we obtain: 

^ <L>=- J f (sAVy)(/) d 3 x — J tp*(x A V F)^jd 3 x . (2.3.11) 

Like in the previous relation, the EIP potential does not contribute to the average 
of the angular momentum. Again if the nonlinear potential has a well behavior at 
infinity the last term is unremarkable on the dynamics of the system. 
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We discuss now the Ehrenfest relation concerning the energy. The energy of a 
canonical system, as we are considering here, is given by E = H where H is the 
Hamiltonian given by Eq. (2.1.31). We can define a Hamiltonian operator H whose 
average value is < H >= H . It is easily verified that: 

H = ~ A + - U Elp (p, j) + - U(p) + V(x, t) . (2.3.12) 
Imp p 

If we compare this expression of H with the expression of the operator O given by 
Eq. (2.3.6) we find: 

H^O, (2.3.13) 

which means that the Hamiltonian operator of a nonlinear canonical system does 
not coincide with the operator O of the r.h.s. of the NLSE whilst, in the case of the 
linear theories, we have O = H = H Q . 

Within the definition E =< H >, we obtain the following relationship: 

which means that when time dependent external potentials are absent, the system 
is conservative being dE/dt — 0. We may conclude that the EIP does not introduce 
dissipative effects. 

We remark that for a noncanonical model, where no Hamiltonian is present, the 
energy of the system is assumed generally as E =< O > where O is the operator of 
the r.h.s. of the corresponding NLSE. Several models with nonlinearities in the r.h.s. 
of the Schrodinger equation, characterized by time independent average values, have 
been developed. For instance, in the Kostin NLSE [87, 88], the operator O is defined 
as O = H + (h r y/2i) [log(ip/ip*)— < log(ip/ip*) >} being a real quantity, the energy 
of the system is defined as E =< O >. In this case the non conservation of < O > 
implies energy dissipation of the system. 

In conclusion, we have shown that the EIP potential (2.1.22) describes a conser- 
vative system. For a free system (V = 0) and when the non-linear potential U(p) 
has a good behavior at infinity, we are able to identify four constants of motion: 

N = J pd 3 x , (2.3.15) 
<P>=-i^J(^*Vilj-ipV^*)d 3 x, (2.3.16) 

< L >= -i ^ J x x (V>* Vip - i> W) d 3 x , (2.3.17) 
E = Jnd 3 x , (2.3.18) 
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representing respectively the energy, the momentum, the angular momentum and 
the number of particles, conserved in virtue of the continuity equation (2.1.15). 



2.3.2 Mathematical proofs 

We deduce here the Ehrenfest relations for the observable N, < x c >, < P >, < 
L >, E given from Eqs. (2.3.2), (2.3.8), (2.3.10), (2.3.11) and (2.3.14). In the 
following we make the hypothesis that the fields vanish steeply at infinity and that 
we can neglect surface terms. 

The first relation, Eq. (2.3.2), is a trivial consequence of the continuity equation, 
while the (2.3.8), can be easily obtained from Eq. (2.1.15) in the following way: 

4- < x c >= 4 [ x^d 3 x = -4 / x V • j d 2 x = 4 / j d 3 x , (2.3.19) 
dt N J dt N J J N J ' 

where we have performed an integration by parts in the last step. 

Let us consider now Eq. (2.3.10) for the component i of the momentum. Using 

Eq. (2.3.4) for < Pi >= i(h/2) f(ip* diip - ip diip*) d 3 x and taking into account the 

following functional derivative: 

- = ihdrt, (2.3.20) 

= Aip + k ^— [if;* Vip - ip ViP*} W 

dip* 2m 2 m 

h 2 

+ k V[ip* Vip-ipVip*} ip + F(p)ip + Vip , (2.3.21) 

Am 

and their conjugate, we obtain: 



j-<Pi> = [{-^-[dtr^ + dnpAip* 
dt J 2m 



h 2 

+ k [ip* Vip - ip Vip*} [Vip* diip - Vip diip*} 

2 m 

h 2 

+ K — V[iP*Vip-ip Vip*} [ip* diip - ip diip*] 
Am 

- [F(p) + V] di p| d 3 x . (2.3.22) 

Integrating by parts twice the first term and one time the third in the right hand 
side of (2.3.22) we obtain: 



^:<p i > = ll-^-d^r^-K^d^rvip-ipvr? 

dt J 2m 8m 
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- [F(p) + V] d iP jd 3 x , (2.3.23) 
and neglecting surface terms we are left with: 

< p. >= J p diVd 2 x + J pdi F(p) d 2 x . (2.3.24) 



d_ 
~di 



Taking into account the relation F = dll/dp, the last integral in (2.3.24) can be 
written as: 

J di{ P F -U)d 2 x , (2.3.25) 

and therefore, if the potential and the field p has a well behavior at infinity, it can 
be ignored. Then Eq. (2.3.24) is equal to Eq. (2.3.10). 

Note that in Eq. (2.3.10) it does not appears the contribution of the nonlinear 
potential U(p). This statement is true also for nonlinear potential dependent on p 
and its spatial derivative U([p\). In fact, in this case we have F = SU/5 p and the 
last term in (2.3.24) is, as before, a surface integral. 

For Eq. (2.3.11) we can follow the same procedure used for the relationship (2.3.10). 
We only give the trace of the procedure. 

Setting < Li >= i{%/2) J e ijk Xj (ip* d k ip — ip d k ip*) d 3 x, by using Eq. (2.3.4) we have: 



< L % > = J e ijk Xji-j^ [dk4>* A-^ + d k ip Aip* 



U 2 

h 2 

+ k — — V [ip* — ip V^*] [ip* d k ip — ip d k ip*\ 
Am 

- [F(p) + V] d kP ^d 3 x . (2.3.26) 

Now we have perform the partial integration, as it was done in Eq. (2.3.22), taking 
into account that the quantity €ij k d k Xj vanishes because the skew-symmetry of e^ k 
and disregarding the surface terms we obtain: 



d ? 

Tt <Lt>= 



J e ijk Xj d k [V + F(p)} p d 2 x . (2.3.27) 



Computing as before the nonlinear quantity F(p) we obtain Eq. (2.3.11). 

Finally, we derive the relation (2.3.14). We begin writing the EIP potential 
as a function of the density p — t/j* ip and of the gradient of the phase S = 
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(ik/2) \og(ip>*/ip>). Equation (2.3.4) becomes: 

;2 ».;.* fc 2 



77 



% z dif)* A , ^ A <9^ 

^ A — - 

2m at 



2m <9t 



d_p , aV 
~di ' 



at ^ EIP ^ dt r dt 
By using the equations of motion of the fields if) and if)*: 



(fx 



i h 



dip 
~dt 



h 2 



2 m 



a^ 



at 



ft 2 



2 m 



A + V + F + W + iW 



A + V + F + W -iW 



if. 



(2.3.28) 

(2.3.29) 
(2.3.30) 



with W, W and F given by (2.1.16), (2.1.20) and (2.1.21) respectively, the Eq. 
(2.3.28) becomes: 



dE 
It 



i^- A if)* A if)- i^— if)* (V + F + W -iW) Aif) 
Am 2 2 m 

i ip* A 2 if) + % — if)* A [(V + F + W + % W) ip] 
Am 2 2 m 



d 3 x 
d 3 x 



+ 



' v ,dU dU^\ dp dU Elp d(VS) 
v dp dp J dt 9(VS) dt 



d 3 x 



dV 

+ < aF > ' 

where we have used the relationship: 



d_ 

dt 



{u + u E 



dU dp dU E1P dp dU E1P d(VS) 



(2.3.31) 



(2.3.32) 



dp dt dp dt d(yS) dt 
In Eq. (2.3.31), integrating by parts and neglecting the surface terms, we obtain: 
dE ih 



dt 



2m 
k 



2 m 



J{V + F + W) {if)* A if) -if) A if)*) d 3 x 
J W (ip* Aip + if)Aif)*) d 3 x 



+ 



, v+ d£ + dU mp 



dp dp J dt 







dS) 


dt 


d(vs)_ 


~dtj 



dV 



(2.3.33) 
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Using Eqs. (2.2.1) and (2.2.2), that we can rewrite in the form: 



d -P = -v 

dt 



'^S \ 2 
m I h 



ds _ h 2 Ayfp (vsy 



dt 2m y/p 2m 
and taking into account the relations: 

dU dU^ 



V -F-W 



dp dp 



F + W 
2 



h 



i h 



d(ys) 



2m 

(^* A^p + ^A^*) = 2p 
Eq. (2.3.33) becomes: 

= j {v + F + w)v{^ P )s 

J(V + F + W) 



m 



dE 
~dt 



x 



h 



m 



+ 



\ 2 yA ' 
m I h 



pW 

d 3 x 



2 r 
+ -JpW 

av 
+ < W > 



r *2 



h Ay/p (ysy 



2m yfp 2m 



V -F-W 



Vp 



d 3 x 



We can immediately obtain: 



dE 
It 



< 



dV 
~dt 



> , 



(2.3.34) 
(2.3.35) 

(2.3.36) 

(2.3.37) 
(2.3.38) 

(2.3.39) 



h 



d 3 x 



(2.3.40) 



(2.3.41) 



that is the relation (2.3.14). We note the absence of contribution of the nonlinear 
potentials C/ E ip(p, j) and U(p) to the average of the energy of the system. We 
conclude by seeing that relation (2.3.41) can be obtained straightforwardly by using 
the following property of the Poisson brackets: 



{/, /} = , 



(2.3.42) 
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valid for every functional / and therefore, we have: 

IT = < aF> • (2 ' 3 ' 43) 

if we take into account that explicit time dependence in the Hamiltonian can occur 
only in the external potential V(t, x). Notwithstanding we have preferred to obtain 
it in a more rigorous fashion because of its importance. 
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Chapter III 

Symmetries and Conservation 
Laws 

The concept of symmetry plays an important role in the search of the solutions 
of dynamical equations. In the case of dynamical systems with infinite degrees of 
freedom described by partial differential equations (PDE), the Liouville integrability 
requires the knowledge of infinite symmetries. In fact many integrable PDE show 
this property like, for instance, the linear Schrodinger equation iip t + ip xx = [66], 
the Korteweg-de Vries (KdV) u t — 6uu x + u xxx = and the modified KdV (mKdV) 
u t - 6u 2 u x + u xxx = [89], the cubic NLSE iip t + ip xx + |^| 2 V = [90, 91], the 
Kaup-Newell equation iip t + i>xx —i(^>^>* ^>) x = [92], the Chen, Lee and Liu equa- 
tion iip t + ip xx — i \ip\ 2 ip x = [26] and others. 

When the evolution equation can be obtained from a variational principle, the pow- 
erful Nother theorem [82] gives the possibility to compute directly the associate 
conserved quantities satisfying appropriate continuity equations. The symmetries 
associated to geometrical transformations are very importants. In fact, as it is 
known, these are related to conserved physical quantities like the total energy, the 
linear and angular momentum of the system and so on. 

In this chapter we study the symmetries described by the connected local Lie groups 
Q depending on r parameters ^(t, x, ip, ip*) which are functions of the independent 
variables t, x and the dependent ones ip and ip*. They are continuous transfor- 
mations mapping solutions of the evolution equations obeying the EIP in other 
solutions. We show that in presence of the potential introduced by EIP some geo- 
metrical symmetries of the linear Schrodinger equation are lost. 
When the full group of symmetry Q is known, we can study the set of solutions which 
are invariant under the action of symmetries belonging to the coset Q/K where /C 
is the set of subgroup of Q in r — 1 parameters. In fact in this case, the original 
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PDE is reduced to an ordinary differential equation (ODE) which can be solved 
by quadrature. We perform this computation in chapter V where soliton solutions 
obeying the EIP are obtained with this procedure. 

At the end of this chapter we will discuss in some detail a particular class of non- 
linear gauged symmetries which can be considered as a generalization of the gauge 
transformation introduced recently by Doebner and Goldin [54]. Because it can be 
applied to a wide class of NLSE, we will study this transformation in a general 
contest. The importance of this transformation is that in order to make real the 
complex nonlinearity that occurs in the NLSE and as a consequence it linearizes the 
continuity equation. In some cases, known in the literature [54, 66, 93, 94, 95] it is 
also a powerful tool to linearize the NLSE at whole. 

3.1 Lie Symmetries 

Let us rewrite the EIP-Schr6dinger equation (2.1.24), after an appropriate rescaling 
of the space-time coordinates, in the form: 

i i> t + i> xx + « (r ^ x - i> r x ) ^ x + \ ir ^ xx - $ r xx ) $ = o , (3.1.1) 

where we use of the stenography notation: ip x = V^, i ) xx = A^. 
We are concerned with the study of the symmetries of the EIP potential and thus, 
we neglect in this chapter both the nonlinear and the external potentials U (p) and 
V. The results obtained can be generalized when these quantities are present. 
Being Eq. (3.1.1) canonical, we can perform the study of the symmetries both at 
the level of the Lagrangian and directly on the evolution equation. The study of 
the symmetries from the PDE is more general, because as we will show, not all the 
symmetries of the PDE are also symmetries of its Lagrangian. Therefore we start by 
studying the symmetries of Eq. (3.1.1) and after in the next section, starting from 
its Lagrangian, we derive, by means of Nother theorem, the appropriate conserved 
quantities (when it occurs). 

To find the Lie symmetries of Eq. (3.1.1) we adopt the geometric method presented 
in the Olver text-book [65] to which we send the reader for a complete treatment. 
We make a summary of the general concepts. 

The first step in the research of the symmetries for a PDE is to look up it as 
an algebraic equation. For this purpose, we introduce the jet-space F = M x 

X XX 

F x F x F x F l which is an extension of the configuration space. Here M 

is the configuration space mapped by the space-time coordinates, F is the space 

mapped by the function tp and tp* when are considered as independent fields. In 
x 

the same fashion F is mapped from the function and tp* and so on. In this 
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jet-space, Eq. (3.1.1) becomes an algebraic equation in the independent variables 
x, t, if), ip*, -0> ip*, V'ttj VC™ an d i>t an d defines an hyper-surface on F. We 
consider now an element g G of a transformation Lie group which acts on the 
element of F. We say that g is a symmetry transformation for Eq. (3.1.1) if a point 
P is mapped on the hyper-surface in a point on itself. In particular, if we consider 
infinitesimal transformations, we must require that those generated by the element 
g lie on the tangent plane on the hyper-surface in P. 

Because Q is a Lie group, its elements around the identity can be written as: 



9 



JeV 



(3.1.2) 



where e is a parameter and v is a vector in the tangent plane on the hyper- surface. 
Its general expression will be: 



v = r d t + £ X ■ d x + + dp + 0* dp 

X „ ~X „ , XX 



+ ■ d^ x + • d^ x + 
where the coefficients: 



d *xx + ^ 



~ £C£C 



a, 



^xx 



JX 




4> 



X 



( V 

<py 



X 



( (f) x \ 
V 4> z / 



(3.1.3) 



XX 



= 




-XX 





I <p xx \ 
\4> zz J 



(3.1.4) 



are three- vectors. All this quantities together with r, and are function of x, t, ip 



rv> . rr» rr» rr» . ry ry 



and ip*. The quantities 0,0 ,0 ,0 and 0* are related to the coefficients 



.x 



(3.1.5) 



£ and r by means of: 

0* = Dj — (d, n ^ - (Dj e) ^ - (Di e) ^ - (d* ^) ^ , 

with i = x, y, z, t where 

Di = <9i + V»i <9</, + ip* dp , (3.1.6) 
is the total derivative, and: 

^ = D,,0 - 2 (D,n ^ - (D.f) V. + (x - y) + (x - z) + (x - t) . (3.1.7) 
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In Eqs. (3.1.5) and (3.1.7) the time derivative of the fields ip and ip* can be eliminated 
by using the evolution equation (3.1.1). Of course, analogous expressions hold for 
* and 11 . At this point we must determine the coefficients r, £ x , and so that 
the vector (3.1.3) lies on the tangent plane to the hyper-surface generated by Eq. 
(3.1. 1) 3 . Let L(ip, V*) denote the PDE: 

L(V, r) = ih + 1> xx +KW1> x -W x )1> x + ^ (r i> xx -i>i>* xx )i> ,(3.i.8) 

g is a symmetry transformation if: 

L(grl>,gP) = 0, (3.1.9) 

whenever ip and ip* are solutions of L(ip, ip*) — 0. The condition for the vector v to 
lie on the tangent plane is given by [65]: 



v 



iipt + ip xx + K(ip*^ x -ip^* x )^ x + ^(^*^ xx ~^V XX )^ 



= . (3.1.10) 



Taking into account the expression of the vector v (3.1.3) we obtain: 

1^ + + 4> xx -^ 2 4> xx + k (2f^-#i) (i> x -K^^x$ x 

+ | (ri'xx -2W XX - 24> x r x ) 0+ \ (Vfe + Wx) = 0, (3.1.11) 

which must be solved using the expression of the quantities 0, 0*, M , * and M . 

By requiring that the coefficient of each monomial in ip, ip* and their derivative 

vanishes separately, we obtain a sequence of derivative equations for the coefficients 
x ~ 

r, £ , 0, that when satisfied, give us the expression of the generators of the full Lie 
symmetries. The computation was developed with MATHEMATICA 2.0 package. 
We report in table the results: 

d xx e+d yy e+d zz e = o, 0»e = -0*e, o t e = 2d x e, j> = ^, 

d xx ^ + d yy ^ + d zz ^ = , d z S* = -d x ?, d x ? = d y ?, 0=-^* ,(3.1.12) 

d xx e+d yy e+d zz e = o, ^e = -9„r, ^e = ^e, 

and all the other derivatives are equal to zero. The general solution of this system can 
be easily obtained. Thus, we conclude that the most general infinitesimal symmetry 



3 In Ref. [65] the generator of the transformation v acting on the jet-space F was called prolon- 
gation of the vector u and denoted by v = pr^u where u is a vector on the tangent plane on M. 
Here we prefer to avoid this terminology for purposes of simplicity. 
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of the EIP-Schr6dinger equation has coefficient functions of the form: 





= A (y — z) + n x + a x , 


(3.1.13) 


? 


= A (z — x) + ji y + a v , 


(3.1.14) 


r 


= A (x — y) + /i z + a 2 , 


(3.1.15) 


T - 


= 2 /it + a* , 


(3.1.16) 


= 


= P1>, 


(3.1.17) 


= 


= -w , 


(3.1.18) 



where A, /i, /3, a 1 are real coefficients. Thus the Lie algebra of infinitesimal symme- 
tries is spanned by the nine vector fields: 



v 2 
v 3 

V A 


= d x ' 

= 8y 

= d z 

= d t j 


> 




translation, 


(3.1.19) 


v 5 
v 6 
v 7 


= yd x 

= xd z 
= zd y 


— xd x 

- zd x 
-yd z t 


> 


rotational, 


(3.1.20) 


v s 


= t d t + 2 (x d x + y d y + z d z ) 


dilation, 


(3.1.21) 


Vg 






U(l) global transformation. 


(3.1.22) 



The groups Qi generated by the Vi are given in the following table. The entries give 
the transformed point exp(eVj) (t, x, ip, ip*) = (t, x, ip, ip*): 



y 1,2,3 


(t, x + e,tp, ip*) , 


(3.1.23) 




(t + e, x, ip, , 


(3.1.24) 


Qb,6,7 


(t,R^x v ^, r) , 


(3.1.25) 


Qs 


(e 2 H, e*x, ij, r), 


(3.1.26) 


G9 


(t, x, e^^e^r) , 


(3.1.27) 



where R l i is the 3x3 relational matrix. Since each group Qi is a symmetry group 
we have that if ip(t, x) is a solution of Eq. (3.1.1), so are the functions: 

= ijj(t, x + e) , (3.1.28) 
^,( 2 ) = ijj(t + e, x) , (3.1.29) 
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V> (3) = i/>(t, R ij 'xj) , (3.1.30) 
^( 4 ) =ij{e 2e t, e e x) , (3.1.31) 
^(5) = e i£ ^(t, x) , (3.1.32) 



where e and e are real quantities. 



3.2 Conserved quantities 

Having obtained the generators of the Lie symmetries and their algebrae, here we 
study the conserved quantities related to these symmetries, starting from the action 
of the system and by means of the Nother theorem. In this section it is convenient 
to restore the standard unity of h and m. As it was said previously, not all the 
symmetries of the PDE (3.1.1) obtained in the last section are also symmetries for 
the action: 

A= I Cd 3 xdt . (3.2.1) 

As a consequence do not all those symmetries are related to conserved quantities. 
To start with, we remember that " Lagrangian" symmetries are those that do not 
change the formal expression of the action. Therefore, they are coordinates and 
fields transformations satisfying the relationship: 



5A = J [5Cd 3 xdt + C5(d 3 xdt)]=0 . (3.2.2 



The Nother theorem states that a current J v exists, given by: 
satisfying the continuity equation: 

d„J v = . (3.2.4) 

Therefore the quantity: 

Q = J J°d 2 x , (3.2.5) 
v 

is time conserved if the spatial component J % vanishes steeply on the boundary dT>. 
Taking into account the result of the last section we begin considering the invariance 
of (3.1.1) under U(l) transformations (3.1.22): 

i//(t, x) = e- ia/n *P{t, x) . (3.2.6) 



41 



The current is: 



cp, 



i h 



2m 



and satisfies the continuity equation: 



(3.2.7) 



(3.2.8) 



The time component of Eq. (3.2.7) is the matter density p, the spatial part is the 
quantum current j with EIP. 

Consider now the space-time translations with generator (3.1. 19): 



t — > t — to , X — > X — Xq , 

the field ip(t, x) transforms as: 

lfj(t, X) — > l/)(t + t , X + Xq) . 

These symmetries imply the conservation law: 

where the energy-momentum tensor T^ u obeys to EIP and is defined as: 



rpOO 
rpOi 

rpiO 



2m 



h 2 



8m 



-i c - (ip* dii/j -ipdii/j*) , 



h 2 



(dill)* doll) + dfll)*) 



rpij _ 



2 m 

h 2 

+k — (r - i) d^*) (r d *p - i> a r) 

Am 



(diip* djip + Oil/) dji/)* 



2m 



h 2 



— K ■ 



Am 



(V>* dii) - ip dii)*) (ip* djip - ip djip* 



-5, 



ti 2 f) 2 

— Ap + K — (iP*Vil)-ipV ip*) 2 
Am 8m 



(3.2.9) 
(3.2.10) 
(3.2.11) 

(3.2.12) 
(3.2.13) 

(3.2.14) 



(3.2.15) 



Equations (3.2.12) and (3.2.13) define the generators of the transformations (3.2.9) 
and are the energy density and the momentum density of the system, respectively. 
From Eq. (3.2.12) we may define the quantity: 



u^ = -K^{rvip-ipvr) 2 , 

8m 



(3.2.16) 
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representing the interaction energy density introduced by EIP. The quantity T 00 
is the sum of the kinetic energy density {h 2 /2 m) |V?/>| 2 and the interaction energy 
density: 

T 00 = ^|V^| 2 + f/ EIP . (3.2.17) 

From (3.2.13) we can see that the momentum density is not influenced by EIP. The 
energy flux T t0 and the momentum flux T y are modified by EIP, as we can see from 
Eqs. (3.2.14) and (3.2.15). Using Eq. (3.2.11) we remark that the quantities: 

E = J T 00 d 3 x , (3.2.18) 

P i = 1 j 'r 0i d 3 x , (3.2.19) 

are constants of motion. By comparing Eqs. (3.2.13) and (3.2.14) we can see that the 
momentum density T 0i does not coincide with the energy flux density T l ° because 
of the non Lorentz invariance of the system. On the contrary we have T JJI = T Jl , 
suggesting spatial rotations invariance of the system. 

In a nonrelativistic theory, like the one we are considering, the energy of the 
system must be a semidefinite positive quantity. This is a true statement for our 
model. In fact, taking into account the expression of the spatial component of 
(3.2.7), which represents the quantum current density, Eq. (3.2.12) of the energy 
density T 00 becomes: 

T oo = ^!_| VV ,|2 + k ™{_3_) 2 . (3.2.20) 
2m 1 ^' 2 \l + K P ) V ; 

This expression appears to be a semidefinite quantity when the parameter k is 
positive. We can transform this expression in the form: 

T 00 = ^|V^(l + M-^(Vp) 2 , (3.2.21) 
2m 8m 

which now results a semidefinite positive quantity when k, is negative, if we remember 
that the quantities 1 + k p is always positive. 

Let us consider the transformations on the coordinate x, produced by the or- 
thogonal matrix R: 

- /.", .'•' • 11^11^ = 6^, (3.2.22) 
ip(t, x) -> ip(t, R' 1 x) . (3.2.23) 
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The action invariance allows us to define the angular momentum density: 



M m = x 3 Tf* - x l T w , 



(3.2.24) 



with ji — 0, . . . , 3; i, j — 1, 2, 3, obeying the continuity equation: 



df 1 M m = . 



(3.2.25) 



We may note from Eqs. (3.2.13) and (3.2.15) that the generators of the transforma- 
tion (3.2.23) 



are constants of motion, not affected by U Elp which, on the contrary, modifies the 
flux densities. 

We are left with the scaling symmetry generated by Eq. (3.1.21). It is easy to see 
that it is not a symmetry for the Lagrangian, independently of the number of the 
dimension D in which the system is immersed. Therefore we are not able to apply 
the Nother theorem to obtain the correspondent conserved quantities. 

We ask now what happens to the other generators of the Schrodinger group. In 
particular, from the list of symmetries obtained in section 3.1, the Galilei and the 
conformal symmetries are missed. 

Let us consider the conformal group transformations: 



with a, (3, 7, 5 arbitrary constants, that can be decomposed in three independent 
transformations: t — > \/t, t — > at and l/t — > t + S. The first is a discrete transfor- 
mation and does not produce constants of motion, the other two, the dilation and 
the special conformal transformation, makes the action of the system not invariant. 
The reason why we loose the conformal symmetry lies in the fact that the parameter 
k has a proper dimension. In order to analyze how the potential U EIP breaks the 
symmetry, consider the case of dilation: 




(3.2.26) 



t = 



rt + s ' 



a 8 — /?7 = 1 , 



(3.2.27) 




(3.2.28) 



In this case, Eq. (3.2.2) becomes: 




(3.2.29) 



where: 



(3.2.30) 
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is the well-known dilation current of Nother. Therefore the potential ?7 EIP prevents 
to write the quantity d^D^ 1 + U EIP in the (3.2.28) as a tetradivergence. Note also 
that the conformal invariance is generally broken when derivative potentials [96], 
like U EIP , are present. Nevertheless, if we consider the transformation: 

[x, t, k, ip(t,x)] — > \x, \ 2 t, \p\~ 2 k, pLt/j(t, x) , (3.2.31) 

with A and p arbitrary constants, Eq. (2.1.25) remains invariant. Remark that in Eq. 
(3.2.31) the transformation with parameter A is the same of the scale transformation 
found in the last section while, the transformation with parameter p is not a Lie 
symmetry. It permits us to reduce the study of the system described by Eq. (2.1.25) 
to consider only the two relevant cases with \k\ = 1. 
We consider now the Galileo transformation on the coordinates: 

*->f, x^x + vt, (3.2.32) 

and we set the field transformation ansatz: 

ip(t, x) -> R(t, x) e ia[t ' x) ip(t, x - v t) , (3.2.33) 

with a(t, x) and R(t, x) arbitrary real functions. The requirement that the ac- 
tion (3.2.1) be invariant under the transformations (3.2.32) and (3.2.33) gives the 
following derivative on the functions a(t, x) and R(t, x): 

. . mv 1 

Va t, a: = — t-t- , 3.2.34 

n 1 + up 

da(t, x) = mv 2 1 35 

dt 2h 1 + Kp ' 

R(t, x) = const , (3.2.36) 

that, in the limit k — > 0, reproduce the factor a(t, x) = —mv ■ x + mv 2 t/2 of the 
linear theory. We can easily realize that in the general case k ^ Eqs. (3.2.34)- 
(3.2.35) are not satisfied by an arbitrary function of p. 

The Galileo broken symmetry due to the potential (3.2.16) is not surprising. In fact 
we know that the generator of the Galilei transformation is given by: 

G=<P>t-mN<x c >, (3.2.37) 

which represent the velocity center of mass of the system. If we derive with respect 
to time Eq. (3.2.37) and remembering that < P > is time independent because is 
a constant of motion, we obtain: 

^=<P>-mN d< j Xc> . (3.2.38) 
dt dt v 1 
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Taking into account the relations (2.3.9) and < P >= m J j d 2 x we have: 

^ = - Km J pj d 2 x , (3.2.39) 

where j is the current without EIP. We see immediately that the presence of EIP 
breaks the Galilei symmetry which is restored in the limit of k — > 0. 

We conclude with a brief discussion of the discrete transformations P and T. 
These are not broken by the presence of the potential (3.2.16), as it can be easily 
shown considering the following relations: 

P : (t, x) -> (t, -x) , ip(t, x) -> i[>(t, -x) , 

j(t, x) - -j(t, -x) , (3.2.40) 

T: (t, x) -> (-t, x) , if>(t, x) -> x) , 

i(t, x) -> -j(-t, a:) • (3.2.41) 

Therefore the expression of the potential (3.2.16) is invariant under the above trans- 
formations. 



3.3 Nonlinear gauge transformations 

In this section we describe a particular class of nonlinear gauge transformations 
allowing us to linearize the continuity equation making real the complex potential 
that appears in the evolution equation (2.1.25). Since the method that we are 
describing can be applied to a wide class of nonlinear Schrodinger equations, let us 
describe it in a general fashion and only at the end of the section we apply it to the 
EIP system. We introduce the method on the more simple 1-dimensional system, 
the generalization to the 3-dimensional case is in progress. 
We introduce the density of Lagrangian 




2 



2m 



dip 



d x 



2 



U([P], M) , (3.3.1) 



which describes a class of one dimensional nonrelativistic and canonical quantum 
systems. In Eq. (3.3.1) the nonlinear real potential U([ip*], is a functional 
of the fields ip and ip*. We will use the hydrodynamic fields p(x, t) and S(x, t) 
[cfr. Eqs. (2.1.5)-(2.1.6)]. The evolution equations of the fields a = ip, ip*, p, S 
can be obtained from the action of the system A — J Cdxdt by using the least 
action principle E a (A) = 0. We make use of the following notation for the spatial 
derivatives: 

d n a 

a " = ^7 ' a = a . (3.3.2) 

ox 11 
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It is immediate to obtain the evolution equation of the field ip that is given by the 
following Schrodinger equation which contains a complex nonlinearity: 



ih ^ = ~^j£ + WM ^ IK < P]) ^ ■ (3-3.3) 

The real W([p], [S]) and the imaginary W([p], [S]) part are given by the following 
expressions: 



W([p], [S]) = j-J U([p], [S])dxdt , (3.3.4) 

mpl [s]) = his I u{[pl [s])dxdt ■ (3 - 3 - 5) 



The evolution equation for the field p is obtained directly from E a (A) = 0, posing 
a = S (now S is the field canonically conjugated to the field — p). We obtain the 
equation: 

§7 + 2& = ^„ (H ,[S]), (3.3.6) 
where the quantum current takes the expression: 

d 



95, 



n+l 



^([P], [5]) 



(3.3.7) 



We remark that Eq. (3.3.6) is the continuity equation and the term in the right 
hand side represents a source for the field p. When the conservation of the number 
of particles N = j pdx is required, the hypothesis that the potential i7([p], [S]) does 
not depend on S but only on its derivative must be introduced, therefore the Eq. 

(3.3.6) takes the form: 

f? + £-°- 

We note that Eq. (3.3.8) can be obtained directly from (3.3.3) and from its complex 
conjugate, performing the standard procedure. We can see that the imaginary part 
W([p], [S]) is responsible for the nonlinearity of the expression of the current 

(3.3.7) . 

Let us introduce the following transformation for the field ip: 

,/>(x, t) - 0(x, t) = U{[ri m) if>(x, t) , (3.3.9) 

which allows us to eliminate the imaginary part of the evolution equation of the 
field ip, which corresponds also to linearize the expression of the current j^. 
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The operator U, generating this transformation, is unitary W = U 1 and is defined 
by: 



wan m: 



ex P 1 tL' 

I ' t n=0 



1 9 
-1)™ / - 



pdx n 



d 



dS, 



n+l 



U([ P ], [S]) 



If we write the field in terms of the hydrodynamic fields p, a: 



4>(x, t) = p 1/2 (x, t) exp 



h 



a(x, t) 



dx\ . (3.3.10) 



(3.3.11) 



and, due to the unitarity of the transformation, the modulo of is equal to the 
modulo of the field ip, while the phase a is given by: 



a = S + m ]T(-1) 



n=0 



1 <9 r 



pdx n 



d 



d S n+ i 



U([ P ], [S]) 



dx . 



(3.3.12) 



By accepting the statement made by Feynman and Hibbs ([97], p. 96): "Indeed all 
measurements of quantum-mechanical systems could be made to reduce eventually 
to position and time measurements" , (see also [54]) the two wave functions ip and 4> 
represent the same physical system and, as a consequence, we can interpret the Eq. 
(3.3.9) as a nonlinear gauge transformation of the function described by Eq. (3.3.3). 
From Eq. (3.3.3) and taking into account the transformation (3.3.9), it is easy to 
obtain the following evolution equation for the field 0: 

tr in +W([p],[S])cf> 



2m dx 2 



— -m 
2 



1 8 n 



d 



dS n 



+i 



U([pl [S]) 







- ™£(-i) 

n=0 



-I 

dt \ 

q an 
n=0 P 



1 d r ' 



p dx n 
d 



d 



OS, 



n+l 



d S n+ i 
U([ P ], [S]) 



U([p], [S]) 

0. 



dx 



(3.3.13) 



Note that the nonlinearity appearing in Eq. (3.3.13) is now real. The continuity 
equation of the system takes the form: 



dp d j, 
dt dx 



(3.3.14) 



where the current has the standard expression of the linear quantum mechanics: 

(3.3.15) 



J^ = —P ■ 



m 



48 



The gauge transformation (3.3.9) and (3.3.10) makes real the complex nonlinearity 
in the evolution equation, and makes non canonical the new dynamical system. 
However, this transformation may be useful to describe the evolution of system by 
means of an equation containing a real nonlinearity. We remark that nonlinear 
transformations have been introduced and used systematically for the first time in 
order to study nonlinear Schrodinger equations as the Doebner-Goldin one in Ref. 
[54]. 

The method here proposed can be found in literature applied to equations describing 
systems of collectively interacting particles. 

We can quote for instance the canonical Doebner-Goldin equation [54, 95] that can 
be obtained from (3.3.1) when the potential U([p], [S]) has the following form: 

U([P\, [S\) = j(PiSi-pS 2 ) . (3.3.16) 

A complex nonlinearity is generated in the evolution equation of the field vp, with 
real and imaginary part given respectively by: 

W([p], [S]) = -mD^- (V) , (3.3.17) 

mp]AS]) = ^J^- (3 - 3 - 18) 
The quantum current takes the form of a Fokker-Planck current: 

J, = —p + D Pl , (3.3.19) 
m 

resulting to be the sum of two terms, the former is a drift current while the latter 
is a Fick current. The generator of the transformation U (3.3.10) takes the form: 

U{[ri M) = exp QmD logp) . (3.3.20) 

This is a particular case of a class of transformations introduced by Doebner and 
Goldin and permits to write the evolution equation: 

th ^ = -2^d^ + 2mD Jn^^ (3 - 3 - 21) 



Equation (3.3.21) was studied in Ref. [98]; after rescaling: a — > y 1 — (2 m D/%) 2 a, 
it reduces to the linear Schrodinger equation. 
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Now we describe the nonlinear gauge transformation applied to the Schrodinger 
equation with EIP. The system that we want transform is gives by Eq. (2.1.25) that 
we rewrite here in 1-dimensional space: 



with: 



and 



Hp, L) = W(p, j.) + iWlfi, L) , (3.3.23) 



= -4 v (if^J ■ (3 ' 3 ' 25) 



Therefore, we introduce the unitary gauge transformation U for the field ip: 

*P(x, t) - <P(x, t) = U([p], [S]) ^(x, t) , (3.3.26) 
which acts on the phase of tp as: 

J- X ^J- X =J- X {l + Kp) - (3 - 3 - 27) 
It is easy to see that U is given by: 

U([p], [S]) =exp ^jp^Ldx^j . (3.3.28) 

The current j., associated to the new field vp, takes now the standard form of the 
linear quantum mechanics: 

j* = ~P, (3-3.29) 

* max 

while the continuity equation is written as: 

% + ^ = . (3.3.30) 
at ox 

The evolution equation for the field is again nonlinear: 
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with 



Hp, L) = Km 



p(l + K P ) 



— K ■ 



4m 



d 2 p 
dx 2 




<p , (3.3.31) 



but now it is a real quantity. Note also that the transformation (3.3.26) does not 
affect the extra nonlinearity F(p) because \ip\ 2 = \<p\ 2 = p. Of course the price that 
we pay to make real the quantity A(p, j) given by (3.3.23) is that the new system, 
described by <p, is noncanonical because of the nonlinearity of the transformation. 
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Chapter IV 

EIP-Gauged Schrodinger Model 



In this chapter we introduce the gauged EIP-Schr6dinger model describing collec- 
tive effects in interacting particles systems with the Lagrangian (2.1.23) coupled in 
a minimal fashion with gauge fields that takes values in the abelian group U(l). 
Many are the applications of NLSEs coupled with a gauge field and found in litera- 
ture. One of the more important example is given by the Ginzburg-Landau theory 
of the superconductivity [99], where the same NLSE of Refs. [39, 100] with nonlin- 
earity oc p is coupled with an Abelian gauge field, the interaction of this one being 
described by the Maxwell Lagrangian [47, 101, 102]. 

In the model studied by us, the interaction of gauge fields are described by the 
more complicated Maxwell-Chern-Simons Lagrangian (MCS). Many are the moti- 
vation to consider the MCS interaction respect to the more easy case when only the 
Maxwell term is present. It is well known that the Maxwell theory could be defined 
in any space-time dimension; the field strength tensor is still the antisymmetric 
one F^, the Maxwell Lagrangian C oc F^F^ and the equations of motion do not 
change their form. The only difference is in the number of independent fields con- 
tained in the theory. Differently, when the system is imbedded in a even space-time 
dimensions, there are two possible expressions for a first order derivative gauge La- 
grangian which are both gauge and Lorentz invariant. The first one is the standard 
Maxwell Lagrangian, while the other, C oc & vaT F^ F (TT , can be shown to be a pure 
divergence and therefore does not give contribution to the motion equations. The 
situation changes drastically when the dynamics of a physical system is developed 
in odd space-time dimensions. In this case the interaction of gauge fields can be 
described either by the Maxwell that also the Chern-Simons (CS) terms. 
As a consequence, because of the presence of the CS term, the model is useful to 
describe physical systems with planar dynamics. At this time the reader might 
wonder if our discussion is merely of academic interest. The answer to the question 
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is negative. In fact, two dimensional physics can occur in our three-dimensional 
world. This is because of the third law of thermodynamics, which states that all the 
degrees of freedom are frozen out in the limit of zero temperature; it is possible to 
strictly confine the electrons to surfaces. Therefore it may happen that in a strongly 
confining potentials, or at sufficiently low temperatures, the excitation energy in one 
direction may be much higher than the average thermal energy of the particles, so 
that these dimensions are effectively frozen out. 

How it was suggested by Wilczek [19] the presence of CS term confers to the sys- 
tem an anyonic behavior [28], that now obeys to a non conventional statistics [27]. 
Therefore field theory in presence of CS coupling can describe phenomenologies in 
which particles or elementary excitations could be anyons. This hope has in fact 
been realized in the case of the fractional Hall effect where the quasi-particles are 
believed to be charged vortices obeying to anyonic statistics [103]. Recent exper- 
iments [104] seem to confirm the existence of fractionally charged excitations and 
hence indirectly of anyons. 

Another topic in condensed matter physics where CS term is believed to be cor- 
rect, is the recently discovered high-T c superconductors, characterized by their two- 
dimensional nature [20, 105]. This hypothesis on the usefulness of the CS term is 
confirmed by the P and T violating symmetries which are observed in this materials 



We remember also that CS term is an alternative method respect to the Proca La- 
grangian for giving mass to the gauge field, without breaking the gauge invariance 
[109, 110, 111, 112, 113]. Moreover, it provides an example of topological field theory 
(for a review see [114]) since even in a curve space-time, the action of the CS term 
has the same form without any additional metric insertions. This fact has important 
consequences because, how we will show, the CS term does not give contribution to 
the energy-momentum tensor of the system. 



Let us introduce the density of Lagrangian of the Schrodinger equation obeying to 
an exclusion-inclusion principle with Maxwell-Chern-Simons interaction: 



[106, 107, 108]. 



4.1 MCS model with EIP 




'gauge i 



(4.1.1) 



where £ mat is obtained from Eq. (2.1.23) after the substitution: 
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and takes the final form: 



2 2m 

h 2 

+k— Up* D^j - ij (Dij)*} 2 - U(ij* ij) -Vip*ip , (4.1.3) 

where we have denoted the spatial component of the covariant derivative with D = 
V — % (e/hc) A, the indices are lower and upper depending on metric tensor in the 
Minkowski space rj^ = diag(l, —1, —1): A 1 = i]^ u A u ; U(p) is an analytic real 
potential function of the field p = \4>\ 2 , V is an external potential, m is the mass 
parameter and k the coupling constant for the EIP potential. The quantity: 

U mp = -k£- W Dij-ij (Dij)*} 2 , (4.1.4) 
o in 

is the EIP potential with minimal coupling [cfr. Eq. (2.1.22)]. 
Of course the model might be study in any spatial dimension D (if D is even the 
CS term is absent), but for the purposes of application that will be developed in 
the following chapters, we focus our attention in the planar case D = 2. Because 
the system is in (2+1) dimensions, the greek indices take the value 0, 1, 2 while 
the latin indices, that assume the value 1, 2, are the spatial ones. We have, rc M = 
(ct, —x), dp = (c^d/dt, V). In (4.1.2) c is the speed of light and e is the 
coupling constant of the Abelian gauge field described by the scalar potential A 
and the vector potential A, with A^ = (A , —A). Moreover, we assume the sum 
convention when the indices are repeated. 
The Lagrangian £ gauge is given by: 

r = -1 f F^ v + - F TtIU A F (4 1 5) 

-•—gauge ^ - 1 [j,v - 1 4 1 iiv i y±.j-.aj 

with F^ v = d fl A u — d u A^. The Levi-Civita tensor e TfJ/I/ , fully antisymmetric, is 
defined as e 012 = 1. The parameters 7 and g in Eq. (4.1.5) give the relative 
weight between the Maxwell interaction and the CS one. The result obtained in 
this chapter, posing g = or 7 = respectively, hold also when the Maxwell term 
or the Chern-Simons one are present alone. Because of the planarity of the system, 
the electric field is a vector with component: 

E i = -d A i -d i A° , (4.1.6) 

while the magnetic field becomes a scalar: 

B = -e ij d i A j . (4.1.7) 
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In terms of the component of the tensor F^ u we have: E l = F 0i , B = —Fi 2 . 
Now we introduce the action of the system: 



A 



J Cd 2 xdt . (4.1.8) 



The motion equations for the matter field ip and for the gauge fields can be 
written as: 

E r A = 0, E All A = 0, (4.1.9) 

where E^ is obtained from (1.3.6) with the substitution ip — > A^. Explicitly, Eq. 

(4.1.9) becomes: 

h 2 n 2 

ihcD iP = D 2 i)-n [fD^-^(M] Dip 

2m 2m 

h 2 

—k - — D Dtp-ip (Dip)*} ij + F(p) i\) + Vip . (4.1.10) 
If we introduce the spatial current: 

j = a + kp) w wr\ , (4.1.H) 

given by J = (1 + k p) Jo, being 

Jo = ~[rrh/;-^ (Dr)} =j -—Ap, (4.1.12) 

2m mc 

the spatial current of the MCS theory without the EIP, it is easy to see that Eq. 

(4.1.10) obeys to the following continuity equation: 

_£ + V-J = 0. (4.1.13) 
at 

Eq. (4.1.10) can be rewritten in terms of the fields p and J in the form: 

h 2 

ihcD 4> = -—D 2 iP + A(p, J)4) + F{p)4) + Vi> , (4.1.14) 
where the complex nonlinearity A(p, J) is given by: 

A(p, r>=K™(^\ -iK±D(-±P-\ . (4.1.15) 
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The motion equations for the fields A^, given from the second of Eqs. (4.1.9), are: 

ld,F^+ 9 -e^F T , = -r , (4.1.16) 
where the covariant current J v is given by: J" = (cp, J). 

Eq. (4.1.16) requires same comments. For g = we recognize the standard Maxwell 
equations with sources. It is well known that this equation does not admit any 
trivial solution also in absence of the matter field J. Differently, when we set 7 = 
the Chern-Simons equations for the gauge fields are obtained. This equations does 
not admit trivial solutions only in presence of matter. In fact how we will show 
in section 4.2, in absence of the Maxwell term, the gauge fields can be expressed 
as nonlinear functions of the field ip, and vanish in absence of the matter field. 
Therefore, when the Maxwell term is present, the gauge fields describe dynamical 
fields with proper degrees of freedom, while in absence of this one, the CS term can 
be see as a constrain which the matter field must obey. Is this constrain that confer 
to the system an anyonic behavior. 
Now, we perform the gauge transformations: 

A^ -> A^ + d^ui , (4.1.17) 

^ _^ e -i(e/ftc)«^ ? (4.1.18) 

where a; is a well-behaved function so that e^d^dyoo = (with = —e Vfl ). The 
Lagrangian (4.1.1) changes as: 

C^C + ^ VT d^{ujF VT ) , (4.1.19) 

with an extra surface term that does not change the motion equations of the fields 
ip and Ap. This property of the system is typical in the presence of CS term and 
continues to be valid also when the EIP interaction is introduced. From Eq. (4.1.16), 
contracting with the differential operator d v , we may obtain again the conservation 
law of the current (4.1.13) in the covariant form: 

d v J v = 0. (4.1.20) 

The time component of Eq. (4.1.16) is the Gauss law: 

7 V -E-gB = ep , (4.1.21) 

which is the expression usually reported in literature, because the EIP does not 
modify the time component of the current. 
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After integration on the whole plane of Eq. (4.1.21), and taking into account that the 
CS term is dominant over the Maxwell term at long distance [51, 115], we obtain 
the important property that every configuration with charge Q = e J p(t, x) d 2 x 
transports also a magnetic flux $ = f B(t, x) d 2 x: 

-g<& = Q. (4.1.22) 

This relation suggests the following interpretation: The system described by the La- 
grangian (4.1.1) can be interpreted as analogous to a system of magnetic monopoles 
in (2+1) dimension obeying to a generalized exclusion-inclusion principle in the con- 
figuration space. 

Finally, the spatial component of (4.1.16) is the Ampere law in (2 + 1) dimensions 
with CS contribution: 

1 VAB-^^-gE* = -J , (4.1.23) 

c at c 

where E* is the dual vector of the electric field with components E* = Ej. 



4.2 Hamiltonian formulation 

The Hamiltonian formulation is crucial for the construction of self-dual solution 
which we have analyzed in the last chapter and for a variety of other purposes. In 
presence of gauge fields the Hamiltonian formulation reserves a special treatment. 
In fact, because of the non vanishing Maxwell term (7 7^ 0), the Lagrangian (4.1.3), 
(4.1.5) is degenerate in the velocity and the system can be described only as a con- 
strained Hamiltonian, the constraint being given by the Gauss law. In its treatment 
we follow the Dirac-Bergmann approach [81]. 

The field canonically conjugated of the field 0, is given by Eq. (1.3.17): 
7r = dC/d(j) where the dot indicates the time derivative. Taking into account the 
expression of the Lagrangian given by Eqs. (4.1.3), (4.1.5) and setting for the 
fields ip, ip*, Ap, we obtain the following expressions: 

= i\l>\ (4-2.1) 

v = -4^' (4 - 2 - 2) 

n M = -i^o + l-eo^. (4.2.3) 
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Eqs. (4.2.1) and (4.2.2) give rise to primary derivative: 

■ h i* 



(4.2.4) 
(4.2.5) 



while Eq. (4.2.3) show that Il vanishes and therefore it is a third primary constraint. 
Performing the Legendre transformation: 



we obtain the Hamiltonian density of the system: 

^mat "^gauge 5 

where the Hamiltonian of the matter field 7Y mat is: 



(4.2.6) 



(4.2.7) 



^mat — 



2 m 



n 2 

Di;\ 2 -K— Dip - ip (Dip)*} 2 + U(ip* tp)+ViP*iP , (4.2.8) 



8m 



while the one of the gauge field 7i 

?2 



gauge 



is: 



H gauge = I (E 2 + 5 2 ) + A (g B - 7 V • E + ep) + 9, (A IT) . 



(4.2.9) 

Let us introduce now the Poisson brackets (cfr. section 1.3) between two functionals: 

5 f{x) 5g(y) 5f(x) 5 g{y)' 



{/(«), <K*)} = J 



+ 



+ 



5 ip(z) 5 ir^(z) Stt^(z) 6ip(z) 
5f(x) Sg(y) 5 f(x) 5 g(y 



d l z 



_5iP*(z)5tt^(z) 5tt^(z)5iP*(z) 
5f(x) 5g(y) 5 f(x) 6 g(y) 



d 2 z 
d 2 z 



(4.2.10) 



5A^(z)5U ll (z) SU^z)SA^z) 
The second class primary derivative (4.2.4) and (4.2.5) satisfy the relation: 

{£i(t, x), 6(t, y)} = -ih6<V(x-y) (4.2.11) 
and can be accommodated by the introduction of the Dirac brackets: 

{/(*), g(y)} D = {f{x), g{y)} + 1 f {f(x), {&(*), g(v)} d 2 z 



J{f(x), g{y)}d 2 z . (4.2.12) 
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Expression (4.2.12) can be simplified if one solves the relation (4.2.4), (4.2.5) for 
and it and treats / and g as functional of ip and ip* only. A straightforward 
calculations give: 



{/(*), 9(v)}d = 



1 



+ 



SJjx) 5g{y) SJjx) 5 g{y) 
Sip(z)Sij*(z) 5^*(z)5ij(z) 
5 f{x) 5g(y) 5 f(x) 6 g(y) 



5A^(z)5U lx (z) 5U lx (z)5A^(z)_ 



d 2 z 



d 2 z .(4.2.13) 



The constraint n = is first class and the requirement of its conservation leads to 
a secondary constraint: 



77 = {ITo, H} = -fV-E-gB-ep = 



(4.2.14) 



where H = JTi d 2 x is the Hamiltonian of the system given by (4.2.7). The secondary 
constraint 77 = does not involve any further constraint since, how is easy to verify, 
{77, H} D = 0. 

The total Hamiltonian is now Ht = H+ J A II d 2 x, where A is a Lagrange multiplier. 
How it is well known, A and n o have no physical meaning, and n o = for all 
the time, while A can take arbitrary values. Accordingly we may drop them out 
from the set of dynamical variables of the models. This can be accomplished by 
discarding the term An o in H?, the only role of which is to let A vary arbitrarily, 
and by treating A as an arbitrary multiplier. As a result, the Hamiltonian of the 
gauge fields H ga , uge becomes: 



n 



7 



gauge 



E + B 2 



(4.2.15) 



and the total Hamiltonian density 7i ma t + T^gauge is given in terms of canonical 
fields ip, ip*, Ai, Hi with i — 1, 2. We note that in TC gaVLg e the coupling constant g, 
introduced by the CS interaction, does not appear. 

So far we have assumed 7 7^ 0. In the pure CS model (7 = 0) the situation is 
somewhat different. The theory has two more second class primary constraint: 



& = ni + ^-,4 2 = o, 

2 c 

£ 4 = n 2 - 1-^ = 0, 

2 c 

resulting from the definition of momenta Hj. As 

x),Ut, y)} = 9 -^{ X -y) , 

c 



(4.2.16) 
(4.2.17) 

(4.2.18) 
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these can be accommodated by modifying the brackets: 



{/(*), g{y)} D = {/(*), g(y)} 



r I [{/(*), &(*)}{&(*), <Kv)} - {/(*), {&(*), <Kv)}] d 2 * 

[{/(as), &(*)} {&(*), <?(?/)} - {/(»), &(*)}{&(*), <?(</)}] d 2 * -(4.2.19) 

By solving Eqs. (4.2.16) and (4.2.17) for IIi and n 2 , this can be trough to the form: 

~5f(x) 5g(y) 5f(x) 5 g{y) 



{/(*), 9(y)} 



+ 



% 

h 

I 



5ip(z)5ilj*(z) 5ip*(z)5ip(z) 

S_f{x) Sg(y) Sf(x) 6 g(y) 
5A 1 (z)8A 2 (z) 8A 2 (z)8A 1 (z) 



d 2 z 
d 2 z 



(4.2.20) 



where / and g are considered as functionals of only ip, ip* , A 1 and A 2 . 
In the pure CS case a further reduction is possible. In fact as we will show in the 
next section, one can explicitly solve the derivative (4.2.16) and (4.2.17) and end 
with ip and ip* as the only canonical variables. In the general case, however, the fact 
that the constraint involves both H and Ai makes such reduction impossible. Using 
the Poisson brackets "naked 1 from the constraint we can evaluate the evolution 
equations of the fields ip and ip*: 



i h 
i h 



~dt 



at 



5H 

5 tfj* ' 
SH 



(4.2.21) 
(4.2.22) 



that are equal to Eq. (4.1.14) and its conjugate. The analogous equations for the 
fields A and IL become: 



d A i 
~~dT 

dt 



5H 
SH 

TX 1 



(4.2.23) 
(4.2.24) 



and correspond respectively to Eq. (4.2.3) and Eq. (4.1.16). 
Let us now define the following quantities: 



N = JiP*^d 2 x , 

< x c >= — J ip* xip d 2 x , 



(4.2.25) 
(4.2.26) 
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< P >= -« I J [i/j* Di/j - i> (Dtp)*] d 2 x , (4.2.27) 

<M>=-i\ J x A[ip* Dip — ip (Dip)*} d 2 x , (4.2.28) 

E = Jnd 2 x , (4.2.29) 

that represent, respectively, the particle number, the position of mass center, the 

linear momentum, the angular momentum and the energy of the system. The time 
evolution of these quantities is given by: 

!? = {'• *>» + i?' (42 - 30) 

for a generic functional f(ip, ip*, A i: II.), where the last term in the right hand 

side takes into account the explicit time dependence of the functional /. After the 
computation we are left with the following Ehrenfest relations: 

4-N = 0, (4.2.31) 

dt y ' 

± Jjd 2 x, (4.2.32) 

(i .. . J fd 2 x- yV Wipd 2 x , (4.2.33) 
^ < M >= J x A / d 2 x - J 4)* (x A W) ip d 2 x , 

^E = Jr^^d 2 x, (4.2.34) 



d 

— < x c > 
dt 

d 

— <P> 



where the quantity: 

f = ep+ ^j A B (4.2.35) 

is the density of the Lorentz force with EIP. We note that Eqs. (4.2.32)-(4.2.34) are 
formally identical to the analogous relations of the theory without the EIP (k = 0). 
The difference is in the expression of the current: J = Jq (1 + k, p). From Eq. 
(4.2.31) we see that the number of particles of the system is a constant of motion 
and from (4.2.34) we see that for external time-independent potential also the energy 
is preserved. Note also that EIP potential does not appear explicitly in the relations 
(4.2.33) and (4.2.34) as shown in section 2.3. The same statement is true for the 
nonlinear potential U(p) when it has a good behavior at infinity. Without the 
external force the linear and angular momentum are constants of motion. In the 
presence of a central force, while the linear momentum is not conserved, we have the 
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conservation of the angular momentum. Finally, it is trivial to note that from Eqs. 
(4.1.22) and (4.2.31) we have the conservation of charge and magnetic flux attached 
to the system. 



4.2.1 Derivation of Ehrenfest relations 

The Ehrenfest relations (4.2.31)-(4.2.34) for the observable (4.2.25)-(4.2.29) can be 
obtained following the computation reported in section 2.3.2 for the analogue quan- 
tity without the gauge coupling. The only quantities that involve more complicate 
steps are Eqs. (4.2.33) and (4.2.34). Therefore we derive in a detailed way only 
these two relations. 

In the following we make the hypothesis that the fields vanish steeply at infinity 
and that we can neglect surface terms. 

The first relation, Eq. (4.2.31), is a trivial consequence of the continuity equation 
(4.1.13), while the (4.2.32) is easily obtained from Eq. (4.1.13) as: 

4- < x c >= 4 [ x^d 2 x = -4 / xV J d 2 x = -J- / J d 2 x , (4.2.36) 
dt N J dt NJ NJ v ; 

where we have performed an integration by parts in the last step. 
Relation (4.2.34) is immediately obtained posing f = E being E the Hamiltonian 
(4.2.29) and taking into account the property of the Poisson brackets {/, /} = 
valid for every functional /. Eq. (4.2.30) becomes: 

that is the relation (4.2.34), if we take into account that the explicit time dependence 
in the Hamiltonian can occur only in the external potential V(t, x). 

Let us consider now Eq. (4.2.33) for the component i of the momentum. Using 
Eqs. (4.2.30), (4.2.13) for / =< P i > given by (4.2.27) and taking into account the 
following functional derivative: 

IL = -i h (D^r, (4.2.38) 
' ' Dj D j i) + k U* D j ip - i{j {D j iP)*] Djtj) 



Sip* 2m J 2m 

h 2 

4m 



+ « - — Dj 



i)* D j ij -ij(D j ij)* ij + E( P )ij + vij , (4.2.39) 
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1L 

SH 



e 

— P , 

c 



= 0, 



(4.2.40) 
(4.2.41) 
(4.2.42) 



we obtain: 
d 



dt 



<P l > 



1L 

~2m 



+ K 



i/j* D j i/j - if) (D j i/j)*] [(Dj^Y Dty - Djip (£>V) 



2m 
ft 2 



+ K ^ — 

Am 



(4.2.43) 



Integrating by parts two times the first term and one time the third in the right 
hand side of (4.2.43) we obtain: 



— < P l >-- 
dt 



/{ 2m \-( D 



Am 



[F(/j) + y] d l p + epF M }d 2 x . 



(4.2.44) 



Starting from the identity: 

(DjDW^y $ + Dty {DjD^y = (D i DjD j ipy $ + ([Dj, £)*] L> j v)* Y> 

+ (^'[Dj, D*]^)* ^ + (DjD j ipy (4.2.45) 

integrating by parts the third term in the right hand side of (4.2.45) and noting 
that: 

[D s , £><] = -li - #A 3 ) , (4.2.46) 



we have: 



(DjDW^y V + Dty (DjD j ^y = ^ (fyA* - OMj) [y>* £> j V> - ^ (D j ip)*] , (4.2.47) 
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and moreover: 

(Dj^y Dty - Djij) (£>>)* - ^'Dj-V + V (£>,■>)*] = -* ^ (djA 1 - &A 3 } p .(4.2.48) 
Inserting relations (4.2.47), (4.2.48) in (4.2.44) we can obtain the relation: 

-4- < P i > = [\ -i^- (djA 1 - &Aj) W LPij) - ij> (LP (1 + up) + epF i0 \ d 2 x 
dt J I 2m c V / L J j 

+ y pd i F(p)d 2 x + J pd i Vd 2 x . (4.2.49) 

By taking into account the relation F = dU '/dp, the last two integrals in (4.2.49) 
can be written as: 

I &(pF-U) d 2 x , (4.2.50) 

and therefore, if the potential and the field p have a good behavior at infinity, it 
can be ignored. Then Eq. (4.2.49) is equals Eq. (4.2.33) if we take into account the 
definition of electric and magnetic fields and of current (4.1.11). 
Eq. (4.2.34) can be obtained following the same iter used to deduce the relation 
(4.2.33), we leave the easy but tedious computation to the reader. 



4.3 Conservation laws 

Let us set, in the Ehrenfest relations (4.2.31)-(4.2.34), the external potential V = 0, 
we can deduce the following conserved quantities: 



d 

It 
d „ 

-rrPtot = 
dt 

d_ 

d~t 

d 



N = , N = J pd 2 x , (4.3.1) 

P tot =< P > + J E A B d 2 x , (4.3.2) 

Mtot = , M tot =< M > + J x A (E A B) d 2 x , (4.3.3) 

^-E = , E = f Hd 2 x , (4.3.4) 

N is the particle number, P tot is the total linear momentum, M to t the total angular 
momentum and finally E is the total energy of the system. 

In this section we want to obtain the above constants of motion (4.3.1)-(4.3.4) by 
means of an approach different from that used in the previous section, i.e. by 
means of fundamental principles and of the Nother theorem. The deduction of these 
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quantities is close to the content of chapter III to obtain the energy-momentum 
tensor for the EIP-Schr6dinger equation, but here, the presence of the CS term 
makes the iter a little more complicate. 

We consider the symmetry related to the invariance of the system over space-time 
translations. If we consider an infinitesimal transformation: 



t — > t — a , x — > x — a , 
it is easy to see that the variations of the fields are given by: 

5 a if) = a» d^, 5 a r = a" d^*, 5 a A u = a" d^A v . 



(4.3.5) 



(4.3.6) 



and the function £ changes for a quantity 5 C = 5^ v d v C, (here there is no sum over 
the repeated index u); where the symbol 8^ v is the Kronecker tensor and the index 
fi selects the variation of C over space translations from time translations. 
Introducing the tensor T^ u , defined as: 



rpl 



v space 



where J7~ s p ace is the Nother current generated from a space translation, and 



rp[)v rrv 

1 — Jth 



time i 



(4.3.7) 



(4.3.8) 



where now J^ nric is the Nother current generated from a time translation. Then, from 
the expression of the density of Lagrangian (4.1.1) and using the Nother formula: 



dC 



dC 



d (d„ip) d (d v ip 
and the definitions (4.3.7) and (4.3.8) we obtain the quantities: 



(4.3.9) 



— \Dif)\ A — k — — [if?* Dif) — ip (Dif))*} 2 + U(p) +eA [> p 



2 m 



8m 



+ 7 ( IF^F^ -F° T d°A T ) - 



9 ( 1 



2 V 2 



T e^A T F^-e 0Xr A x d*A 1 



(4.3.10) 
(4.3.11) 



rpiO 



2m 



(D l if))* d°if) + D l if)d"if) 



|0„/,* 



— K ■ 



4m 



- -fF iT d°A T - ^e iXr A x d°A T , 



if)* D> - if; (£>>)*] (ip* d°ip - if) <9>*) 

(4.3.12) 
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2m 



7 F lT 9 J A r - f e lAr A A &> A T + 5 l Uic'^ [r D i; - V [D^f] - ^- \D^\ 



4m 



2m 



+ « A- [V* £>V - ^ W)f - - J ^ ^ + f ^ ^ F, u 



(4.3.13) 



In (4.3.10)-(4.3.13) the gauge invariance of the tensor T^ v does not appear explicitly. 
We consider now the new tensor T ^ given by: 



(4.3.14) 



where the quantity: 



fi»> = A V [ 1 8 T F TfM + 9 - e" Ar F\ T --J» 

2 c 



(4.3.15) 



vanishes because of the motion equations (4.1.16). We observe that the tensor 

T^ v is gauge invariant. Finally, if we remember that the expression of the energy- 
momentum tensor is defined modulo a quantity d T X T ^ U where X T ^ V is a third rank 
tensor which is antisimetric in the two indices r and \i given by: 



X T ^ = - 7 F T » A v + - e^ A A\ A v , 



(4.3.16) 



it is easy to verify that the tensor T^ u defined as: 

= f ^ - 7 8 T (F T ^ A") + - e T » x 8 T (A v A x ) 
is gauge invariant and obeys to the continuity equation: 



(4.3.17) 



<9 M T^ v = d^T^ = . 



As a consequence of (4.3.17) the motion constants of the system: 



(4.3.18) 



(4.3.19) 



are not modified by the substitutions (4.3.14) and (4.3.17). The final expression of 
the components of the energy-momentum tensor T^ v are: 



T 



oo 



±-\D4>\ 2 - « f- [P D i> - i> (Dm 2 + U{p) 
2 m 8m 
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+ 



2 p. . mi 

4 r% i r 



T 0i = ic- 



iP* D> - ip (D V)* 



7 



(4.3.20) 
(4.3.21) 



T 



t0 



2m 



4m 



^* £>ty - ^ (-DV)*] #V - V> (-^V)*] - 7 F ° j ^4.3.22) 



2 m 
ft 2 



— K • 



4m 



ij>* Dty - ip (Dty)* ip* D j ^j -ij(D j ip) 



4 m 8m 



+ U(p)-p 



-o*M}->yF»Ft x -ls*F>"Fr 



(4.3.23) 



We note that in the above expression of T^ v the CS contribution does not appear 
explicitly. This is due to the topological nature of the CS interaction. How we 
have noted in the introduction of this chapter, independently of the geometry of 
the variety in which the system is imbedded, the expression of the CS Lagrangian 
appears the same without the necessity to introduce the metric tensor g^. In fact 
the quantity e^ u \ is a good tensor. As a consequence, if we take into account that 
the tensor can be also obtained by a functional variation of the Lagrangian 
density with respect to the metric we obtain immediately the statement that the 
contributes of CS to the energy-momentum tensor are null due to the independence 
of its from g^ u . Notwithstanding the presence of the CS term modifies the form of 
the gauge fields because they must satisfies Eq. (4.1.16). 

The tensor T^ v is symmetric only in the spatial indices T lj = T j \ because the 
theory is not Lorentz invariant but only rotation invariant and satisfies the continuity 
equation: 

d li T IMV = 0. (4.3.24) 
Of course, the two following quantities are conserved: 



E 



/' 



T 00 d 2 x, 



tot 



T 0l d 2 x. 



(4.3.25) 



In section 3.2 we have shown that the energy density T 00 is a semidefinite positive 
quantity. This properties hold also in presence of the gauge field. In fact, in the 
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case k > by using Eqs. (4.1.11), (4.3.20) and the expression of the fields E and B 
as functions of the potential A^, the quantity T 00 assume the form 

r* = ^\D^ + K | (^f + &(,) + 2 (tf + . (4.3.26) 

Therefore, if U(p) = 0, T 00 is a semidefinite positive quantity. In the case k < 0, 
after rewriting the quantity T 00 as: 

T 00 = £-\Drl>\ 2 (1 + «P) - « ^ (Vp) 2 + U(p) + ]-(E 2 + B 2 ) , (4.3.27) 
2m 8 m 2 

we can obtain the same conclusion, if we take into account that 1 + k p > 0. 

We consider now the invariance of the system under spatial rotations defined by the 

t r ansf ormat ions : 

x i -> tt l 3 x j , Q j k = S\ . (4.3.28) 

The angular momentum density M° and the relative flux M % form the vector M M = 
(M°, iW), that is obtained from the tensor T^ v by the relation: 

Af = eij T« . (4.3.29) 

From the symmetry = T 3 ' 1 we infer the continuity equation for M M : 

d M M" = . (4.3.30) 

The total angular momentum: 

Mtot = ~ c J M°d 2 x , (4.3.31) 

is time conserved and, like Ptot> is n °t modified by the presence of EIP potential. 
The quantity (4.3.25) and (4.3.31) are the generators of the space-time transla- 
tion and rotations supplemented by local gauge transformations when performed on 
gauge not invariant fields. In fact if we consider for example the quantity: 

g = ei aiFi , (4.3.32) 

which generate a space translation with parameter aj, the linear part of the variations 
of the fields can be performed trough the "naked" Poisson brackets (4.2.13): 

ip(t, x) -> ip a (t, x) = gip(t, x) =ip(t, x) + CLi{lp(t, x), P l } D 

= ip(t, x) + a i d i ^(t, x) + ieaiu(t, x)ij){t, x) , (4.3.33) 
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i>*(t,x)^i>* a (t,x) = g*(;*(t,x)=i>*(t,x) + a i {4;*(t,x),P i } D 

= i/j*(t,x) + a i d i ijj*{t,x) + ieaiuj(t,x)ip*{t,x) (4.3.34) 



A(t, x) A(t, x) = gA(t,x) = A(t,x) + a i {A(t,x),P i } D 

= A(t,x) + a i d i A(t,x) + iea i Lj(t,x)A(t,x), (4.3.35) 

Tl(t, x) Il(t, x) = gU{t,x) = U{t,x)+a i {U(t,x),P i } D 

= U(t,x) + aid i n(t,x)+iea i uj(t,x)n(t,x) (4.3.36) 

with uj the gauge parameter. The fact that the quantities (4.3.25) and (4.3.31) 
generate translations supplied by gauge transformations is due to structure of the 
Poisson brackets (4.2.13) which are obtained after taking into account the constraint 
(4.2.14) which is the generator of the gauge transformation (see below). However, 
if we restrict ourselves to gauge-invariant fields, e.g. p(t, x), E(t, x), and B(t, x), 
Eqs. (4.3.25), (4.3.31) generates pure translations: 

{P(t,x),p(t,x)} = Vp(t,x), (4.3.37) 
{P(t, x), F^(t, x)} = VF^a;). (4.3.38) 

It is known that the introduction of the CS interaction confers to the matter 
field a non conventional statistical behavior [51, 112]. This statistical behavior 
is the result of the Aharonov-Bohm effect [116], because, in presence of the CS 
term, each charge field carries forward, intrinsically, also point magnetic vortices. 
This statement is true also in the presence of the Maxwell term, because only the 
asymptotic behavior is invoked, which, because of the high-derivative terms of the 
Maxwell coupling, is dominated by the CS term. This anomalous situation in not 
modified by the presence of EIP because its does not modify the expression of < P > 
and < M >. In fact, posing for simplicity 7 = 0, the second relation of (4.3.25) and 
(4.3.31) becomes: 



<P>=- 1 — J{^j*V^-^V^*)d 2 x-^ J Apd 2 x , (4.3.39) 
< M>= -y J x A (V>* W - i> W) d 2 x - - J x A Apd 2 x . (4.3.40) 
Eq. (4.1.21) can be solved without ambiguity in the Coulomb gauge (V ■ A = 0): 

A(x) = — / V • arctan ( X2 ~ y A p ( y ) d 2 y , (4.3.41) 

2 7T o J \x^ -y^ 
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and inserting it into Eqs. (4.3.39) and (4.3.40), after integration, we obtain: 



<P>=- 1 — J(^j*Vij-ijVij*)d 2 x , (4.3.42) 
<M>=-y Jx A (fV^-^Vf)A- . (4.3.43) 

From Eq. (4.3.43) we can see that the angular momentum of the field ip is the 
sum of two terms: the former represents the orbital angular momentum, while the 
last, the spin, is responsible of the anomalous behavior of the system. This result, 
already known, is not too surprising: its origin is in the Gauss law (4.1.21) that is 
not modified by the presence of EIP. As a consequence of Eq. (4.3.43) it follows 
that EIP does not change the anyonic behavior of the system and so that the spin- 
statistics relation of the anyonic systems holds. Therefore the Lagrangian density 
(4.1.3), (4.1.5) describe a model of anyonic particle with kinetic obeying to the EIP. 
This one generate an exclusion-inclusion effect in the configuration space wilts the 
CS term is responsible of the not conventional quantum statistics. 

Finally, we discuss briefly the gauge transformations. In section 4.1 we have seen 
that the Lagrangian (4.1.3), (4.1.5) is invariant up to a divergence over global gauge 
U(l) transformations: 

ip(t,x) -> e- i{e/hc)ul ilj(t, x) , (4.3.44) 
A^t, x) -> Ap(t, x) , (4.3.45) 

with 8C — 0. The Nother current is now: 

^ = |cp, -^(l + K p)[^*rty-^(rty*)]} , (4.3.46) 

which satisfies the continuity equation (4.1.20). 

We consider now the local gauge U(l) transformation: 

i/>(t,x) -> e- i{ - elhc)w ^ x) ij}{t, x) , (4.3.47) 
A„(t, x) -> A^(t,x)+d^u;(t, x) . (4.3.48) 

It is well known, that the Nother theorem, in the case of symmetries related to 
continuous transformations, does not give conserved quantities but from it we ob- 
tain the same identities. In fact, the gauge charge associated with this continuous 
symmetry vanishes identically. 

In the case of the transformations (4.1.17) and (4.1.18) we have (using the notations 
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of section 3.2): 

5i/j = -^ipu , (4.3.49) 
Tic 

5A^ = d^uj, (4.3.50) 
/" = 9 -^ VT F VT u . (4.3.51) 

The Nother vector = {J°, J % ) can be calculated by using (3.2.3), we obtain: 

J° = (ep-iV -E + gB)u , (4.3.52) 
J l ={^-J' l -1 dpF* - | e^F^ u . (4.3.53) 

By means of Eq. (3.2.4) the charge Q w , given by Q w — J Jq d 2 x, takes the expression: 

Q w = J '( 7 V • E - gB -ep) cud 2 x . (4.3.54) 

We can easily see, after using the Gauss theorem (4.1.21), that is identically zero 
. This means that the Gauss condition (4.1.21) is satisfied at each time, during the 
motion. 
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Chapter V 



Canonical Systems 
Obeying to EIP: Solitons 

In the previous chapters we have considered the definition of collectively interacting 
particles obeying to EIP and we have studied the main properties of a particles 
system. 

In this chapter and in the following one, we study a special class of states of these sys- 
tems called solitons. Solitons are special solutions of nonlinear evolution equations, 
not obtainable within perturbative methods, that preserve their shape during the 
propagation. Soliton solutions appear in many topics in physics. For instance, the 
cubic Schrodinger equation solutions describe propagation of deep water waves and 
of modulated ion-acoustic waves in plasma, three-dimensional diffractive patterns of 
a laser beam [60] and more importantly describe the recently observed Bose-Einstein 
condensation in rarefied vapors of metal 7 Li, 23 Na 87 Rb [42, 117, 118, 119]. 
Due to their not spreading property and to their particle-like behavior, solitons can 
be used in field theory in order to describe elementary particles. Rigorously speak- 
ing, solitons are wave packets which preserve their shape during their propagation 
and after collision with other solitons. Usually it is required that solitons preserve 
asymptotically their speed after collisions, but this is hard to met with the interpre- 
tation of solitons as particles. 

Typically, in literature, soliton term is misused to denote solitary waves. These are 
also wave packets preserving their shapes but no assumptions are made about their 
behavior after collisions. Thus, a single soliton is a solitary wave but solitary waves 
are not necessarily solitons. In order to see if a solution is a "genuine" soliton it is 
necessary to follow it after collision with another soliton requiring the knowledge of 
multi-solitons states. This solutions generally are very hard to be obtained. 
In the following, we use the word soliton but all our solutions are simply solitary 
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waves. 



5.1 Solitons 

We study a particular class of solutions of Eq. (2.1.25) in the free case, i.e. when 
V = 0. In this situation, we can consider the motion of the mass center on a straight 
line with uniform velocity (for a discussion of solitons in an external potential see 
for example Ref. [120]). In addition to the condition V = 0, we limit our attention 
to the one dimensional case when the EIP holds. The extension to highest dimen- 
sions is not straightforward. In particular the easy case of waves with amplitude 
modulated only in one space-dimension is non physical situation because they carry 
infinite energy. 

In chapter III we have studied the symmetry Lie group of the system obeying to the 
EIP. The result was that Eq. (2.1.25) is invariant over roto-translation transforma- 
tions, scaling transformations and the global unitary U(l) group transformations. 
We have a nine parameter full symmetry group. We can use this information to 
search special solutions obeying to EIP. In fact, we can require that the solution 
should be invariant under a proper selected subgroup of the full symmetry group 
so that the solution itself can be obtained solving an ordinary differential equation 
rather than a more complicate partial differential equation. 

In this section we focus our attention on solutions that are constants in between the 
y and z spatial direction and moreover are invariant over the x and t translation: 

x^x±ue, t—*t + e, (5.1.1) 

with e and u constants. The global invariant of this transformation is the quantity 
x^fut. Therefore we require that the field ip depends only on the time t and on the 
coordinate of the soliton mass center £ = x =F u t where u now has the meaning of 
velocity of the soliton (as usual, the sign minus stands for a soliton moving from the 
left to the right side of the x axis, while the sign plus stands for an antisoliton moving 
in the opposite versus). Thus the wave function becomes ip(x, t) = ■?/>(£, t) where we 
include a possible explicit time dependence in the phase. (Rigorously speaking, this 
solution is not invariant over the transformation (5.1.1) because the time dependence 
on S. Notwithstanding, this assumption does not modify the general procedure to 
find soliton solutions, allowing to obtain a more general expression). In the following 
we use the method described in ref. [121], valid for the NLSEs that are most 
frequently encountered in physical problems. We assume that, for £ — > ±oo, the 
particle density — > so that d£ = N, where N represents the collective 
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particle number contained in the soliton. Moreover, the phase S is written as: 

S = s(0-et , (5.1.2) 

and the field ip assumes the form: 

^ = p(0 1/2 exp{l[ S (0- e t]} • (5.1.3) 

It is now easy to verify that the Hamilton- J acobi equations (2.2.1) and the continuity 
equation (2.2.2) describing the solitonic state can be reduced to the following system 
of coupling equations: 



±u 



ds 1 + 2 up ( ds s 



2 m 



+ U q (p) + F(p)-e, 



dp Id 
<9£ m <9£ 



3s 



(5.1.4) 
(5.1.5) 



where we have indicated with U q (p) the one dimensional quantum potential in the 
£ coordinate: 

h 2 l d 2 ^p 



Uq{p) ~ 2m de ' 

The quantum velocity v q (£) = m _1 <9s(£)/<9£ must be finite when £ 
this condition Eq. (5.1.5) can be integrated a first time, obtaining: 

d s mu 
<9£ 1 + Kp 7 

and after second integration with the condition £(0) = 0, we have: 



(5.1.6) 
±oo; with 

(5.1.7) 



S ^ = ±mu I TV^pW) ■ 



(5.1.8) 



Eq. (5.1.2) and (5.1.8) allow us to calculate the phase S(£) provided that the 
quantity is known. 

To evaluate the density we note that, if we take into account Eq. (5.1.7), 
Eq. (5.1.4) reduces to the following second order differential equation: 



2 d 2 p 



' 1 dp" 



+ 



{2mu/h) 2 8m 

(1 + K P ) 2 ~lf 



F(P) + 



8me 



o . 



(5.1.9) 
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Before solving this equation, we observe that e can be written in the form: 

£= ^ + F(p) _l m „ 2 _L_. (5.1.10) 

Equation (5.1.10) has an immediate physical interpretation: the quantum potential 
causes the spreading of the ordinary Schrodinger wave packet; this spreading is 
compensated by the nonlinearity F(p) and by the EIP contribution — (mu 2 /2) (1 + 
up)' 2 . Therefore it is possible to build up a non-spreading solitary wave. 
After the introduction of the function 



so that 



(5.1.12) 



dy _ 2 d (I dp\ 
dp p di \p d£) 

Eq. (5.1.9) reduces to a first order linear differential equation: 

dy y (2mu/h) 2 8mF(p) 8mel 

h — 7T, \2 ~ T2 h_ T^ — = °> (5.1.13) 

dp p p(l + Kp) 2 h p h p 

that can be easily integrated, giving: 

A 8me (2mu/h) 2 8m U(p) /r . , , iN 

2/ P = p- " i , + "T2 , 5.1.14 

p ft 1 + k p h p 

where the integration constant is A— 1/k in order to obtain the right limit for k — > 0. 
By comparing Eq. (5.1.11) to Eq. (5.1.14), we obtain: 

(dp\ 2 . 8me , (2mu\ 2 p 2 8m ~. . . „ . 

m =Ap -^r" -{— ) Trr f + W u(f) - (5 ' L15) 

The evaluation of the soliton shape is thus reduced to the solution of the first order 
ordinary differential equation (5.1.15). By introducing the dimensionless variables: 

n— \k\ p , (5.1.16) 

and 

2 mu , , . _ „ „ N 

X = ^ (x^f ut) , (5.1.17) 
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Eq. (5.1.15) takes the form: 



- — = a n + P n 2 + 7 nU(n) 



an 



(5.1.18) 



where a = A\k\ (h/2mu) 2 and /3 = —2e/mu 2 are the new integration constants, 
7 = 2 I as: I /mu 2 and £/(n) = U(p)- The parameter <r = k/\k\ assumes the value +1 
when the inclusion principle holds (k > 0, n > 0). Accordingly, for the exclusion 
principle (/t<0, 0<n<l)we have a — — 1. Here we have made use of the scaling 
properties (3.2.28) which permit us to take into account only the two special case 
k = ±1. We note, finally, that while solving the Eq. (5.1.18), we have to take into 
account the two arbitrary constants a and (5 that define a family of solutions. 

To determine the soliton solutions we must search the solutions of the first order 
differential equation (5.1.18), varying the arbitrary constants a and j3, while the 
sign a — ±1 is kept fixed. Eq. (5.1.18) after integration gives: 



±X 



n , 



a n + (a a + (3 - 1) n 2 + a (3 n 3 + 7 n (1 + an) U" 



1 + an 

From Eq. (2.3.16), we have for the soliton case: 

< P >= ±Mu , 

where M is defined by: 

+00 

P 



M 



m 



dx . 



1 + KP 

By using Eq. (2.3.18) we can write the soliton energy as: 



dn . (5.1.19) 

(5.1.20) 
(5.1.21) 



E = 



<P 2 > 
2m 



with 



+ — mu 



<p 2 >=h 



p 



1 + Kp 

dip 2 



+ U(p) 



dx 



dx 



dx 



To evaluate < P 2 > /2m we take into account (5.1.3): 

32 ^ *2 - , / . X 2 / , x 2 



<P 2 > 

2m 



= t/KS dx+ J&) pdx 



4 



1 (dp' 
p [dx 



dx + (muY 



(1 + M 2 



dx 



(5.1.22) 



(5.1.23) 



(5.1.24) 
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using Eq. (5.1.9) we obtain: 

<P 2 > 



2 m 



= mu 



[l + np] 2 



dx 



p d ^A dx + eN _ (5.1.25) 
dp 



Considering that u =< P > /M and Eq. (5.1.25), the energy (5.1.22) satisfies the 
following soliton energy-momentum dispersion relation: 



E 



<P>< 



+ 0O 



+ 



/ 



2M 

U(p) - p 



(1 + K P ) 2 

dU{p) 



dx 



P 



1 + up 



dx 



dp 



dx + eN . 



(5.1.26) 



If we chose the constant e, appearing in the phase of ip, as: 



e N 



+oo 

/ 



dx mu 2 

2 



(1 + np)' 



dx 



(5.1.27) 



the energy-momentum dispersion relation assumes the expression: 

E - <P > 2 



2M 



(5.1.28) 



which is related to a free particle of mass M, traveling with momentum < P > = 
±Mu. 

Let us study Eq. (5.1.18) defining the shape of the solitons. We discuss the case 
U(n) = and f3 = 1. As an example, we consider the case where the EIP is reduced 
to the inclusion principle (boson case) and the particle interaction is attractive 
(er = 1). When a = we can obtain the expression of the soliton explicitly. In fact, 
Eq. (5.1.19) now becomes: 



'1 + n 



rtf' 



dn , 



±x = / 

and after integration we obtain the following implicit solution: 

arcothi 



(5.1.29) 



1 + n 



n 



'1 + n ,1 
n 2 



(5.1.30) 
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The value n(x), at the origin n(0) = n , is the solution of the transcendent equation: 



tanh 



1 + n Q 

n 



n. 



1 + n 



(5.1.31) 



with n = 2.27671. We note that for x — ¥ ±oo, the asymptotic form is n ~ 4/x 2 - 
A soliton with this behavior for x ~^ i 00 was taken into account recently in Ref. 
[122]. We note also that at the origin the soliton has an angular point because its 
first kinetic is discontinuous: 



dn 

Ax, 



±1 



n + 1 



The wave function of the soliton is: 

i>(x, t) = ex P |-» 

where 

e = —mu 2 



n(x) 



\ n(x) 



'l + n e t 
n Q h 



(5.1.32) 



;i + n ) 2 ' 



(5.1.33) 



obtained using Eq. (5.1.27), so that the soliton has a particle-like behavior. The 
expression (5.1.32) of the soliton can be used to calculate explicitly the quantities 
N, < P > and E. Because of the Eqs. (5.1.27) and (5.1.33) the number N is: 



N = 



h 



Kmu 



[l+n ) 



1 + n 



(5.1.34) 



On the other hand if we define the mass M as in Eq. (5.1.21): 

2 



M 



1 + n 



mN , 



(5.1.35) 



the momentum is given by < P >= Mu and the energy is E =< P > 2 /2 M. 



5.2 Effective potential 

In this section we prove that the unitary nonlinear transformation (3.3.9) t) — > 

<f>(£, t)) studied in section 3.3, when applied on the solitonic states ■?/>(£, t) exists and 
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that the new states t) are solutions of a Schrodinger equation with an algebraic 
real nonlinearity. 

Let us consider the unitary transformation: 



where U(£) is given by: 



W(0=exp{^[±m«e-s(0]} • 



(5.2.36) 



(5.2.37) 



The new wave function t) has the same amplitude of the wave function t) 
but a different phase: 



<KZ, t)=p(0 1/2 exp{^(±m M e-et)} 



The unitary transformation can be rewritten as: 



V>(£, t) = exp 



where 



r(f) = ±mu 



(«-/r 



After deriving (5.2.39) we obtain the following relations: 

-ir/h 



~dt 

w 



_ i_dv_ 

~dt ~ H~dt 

i d 2 T 1 



1 OT 



2dVd6 d 2 - 



h 2 \d£ 



hd£d£ di 2 



By considering (5.1.2), (5.1.7) and (5.2.40) we obtain: 

dS 1 8T 



and then the current j, given by (2.1.14) becomes: 

dl 



i , dr 

j = ~ZZ l 1 + K P) 177 • 



(5.2.38) 



(5.2.39) 



(5.2.40) 



-iT/h 



(5.2.41) 



.(5.2.42) 



(5.2.43) 



Km 



(5.2.44) 
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In the case of soliton states the Schrodinger equation (2.1.25) is: 

d^j h 2 d 2 ^ 

+ ^-fr^-V^-^Mf (t^V (5-2-45) 
By using (5.2.41), (5.2.42) and (5.2.44) we may write Eq (5.2.45) in the form: 

9r 90 ^^L^A Hl^]La 2 + k p ( dr \ 2 a 
— < p + lTl — - — + — — + ^ 

ih or d<j) in d far \ 1 _ N ± . n Aets 

m d£ d£ Imp <9£ \0f J 
We use now the relation dT/dt = ^fudT/dt; with: 

d JL = ±mu Jl^ (5.2.47) 
(easily derivable from (5.1.7) and (5.2.44)), Eq. (5.2.46) and: 

<H = o — ~xT ' (5.2.48) 



2mp d£ <9£ m d£ d£ 2m \ <9£ <9£ J d£ p 
If we take into account that j<p = ±u p we arrive to the following NLSE 



! R af = "as 8? + 2 m * K " (7T77F * + F( " ) * ' (5 ' 2 ' 49) 

Let us now introduce the variable x and define F e g(p): 

F eS (p) = F(p) + \ mu ^P ^Xlly ' ( 5 - 2 - 50 ) 
then Eq. (5.2.49) can be rewritten as: 

+ d<f> h 2 d 2 <p . . . 

By means of F e s(p) = dU e s(p)/d p we can introduce the potential U e s(p) which is 
given by: 

U cS (p) = U(p) + \mu 2 K— P - . (5.2.52) 

2 1 + k p 
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We remark that the term (nmu 2 /2) p 2 /(l + up) originates from the £/eip(p, j) and 
represent the EIP effect on the shape of the soliton. We remember also that the 
transformation (5.2.36) is noncanonical. In the case of soliton solution (and only for 
this special solution) we are able to write a new evolution equation derivable by a 
variational principle, starting from the density Lagrangian: 

£cfr = Co — U e s , (5.2.53) 

where £ is the density of Lagrangian of the linear Schrodinger equation. Because 
of the noncanonicity of (5.2.36) U e g is not obtainable directly from U^ip- 

Let us consider the transformation recently introduced by Doebner and Goldin 
[54, 95]: 



ip(t, x) — > <fi(t, x) = y p(t, x) exp 

Transformation (5.2.54) defines a class of non-linear gauge transformations, varying 
the parameters "f(t), X(t) and 0(t, x) and has the important property of making lin- 
ear a particular sub-family of equations belonging to the Doebner-Goldin equation 
family. By comparing Eq. (5.2.37) to Eq. (5.2.54) we can note that the transforma- 
tion introduced in this work is a particular case of the more general transformation 
introduced by Doebner and Goldin. Transformation (5.2.37) is limited to the soli- 
ton states without linearizing the Schrodinger equation describing these states. It 
reduces the complex nonlinearity of the Schrodinger equation to another real one. 
By following the same procedure of the previous section, we can determine the shape 
of the solitonic solutions of Eq. (5.2.51), that are the same states of Eq. (2.1.25) 
modulo the phase S(£, t). Of course, standing on the unitarity of Eq. (5.2.37), the 
equation for the shape of the soliton is the same of Eq. (5.1.18). The problem of 
searching the solitonic solutions of Eq. (2.1.25) with derivative complex nonlinear- 
ities due to the EIP, is reduced now to the search of the solitonic solutions of a 
Schrodinger equation with analytic real nonlinearity. 
We consider now the nonlinear potential: 

U(p) = U Q (p)-\ mu 2 K-^- , (5.2.55) 

I 1 + K p 

where Uo(p) is again an analytic real arbitrary potential in p and the second term is 
selected with the scope of eliminating the effect of the EIP. Eq. (5.1.18) now takes 
the form: 

( dn\ 2e 2 2 \k\ - 

— =an -n 2 + ^rnU (n) . 5.2.56 

\axj mu z mu z 



m {og p ^ x) + J± S (t, x) + e(t, x] 



(5.2.54) 
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Equation (5.2.56) is identical to the equation of the solitonic shape that we may 
find in literature if we take a NLSE with the analytic non-linear potential Uo(p) 
[121]. Then Eq. (5.2.56) allows us to use the soliton solutions of NLSEs available 
in literature. 



5.3 Applications 

As a first application of the results obtained in the previous section, we derive 
the nonlinear potential U(p) which, when is present simultaneously with the EIP 
potential U EIP (p, j), permits the formulation of a soliton with shape given by: p(£) oc 
[cosh(o£)]~ 2 . We start considering the non-linearity [39, 123]: 

U (p) = ^p 2 . (5.3.57) 

The Schrodinger equation with this non-linearity has been recently used to study 
the Bose-Einstein condensation (BEC) [41, 42, 43, 44, 45]. In the Gross-Pitayevski 
equation the parameter p is given by: 

4nh 2 Na 

p = , 5.3.58 

m 

where N is the number of atoms in the condensate, m their mass and a is the s-wave 
triplet scattering length. Its value is assumed to range in the interval 85 a < a < 
140 do, oo being the Bohr radius [124]. 

Set a' — and p > 0, Eq. (5.2.56) with the potential (5.3.57) is easily integrable 
obtaining: 

p(0 = ^ [cosher 2 , (5.3.59) 



b = , (5.3.60) 



where 6 is a dimensionless constant defined as: 

pmN 

~2h 2 

and the normalization N = J l^l 2 ^, that fixes the parameter e = —p 2 mN 2 /8h 2 : 

has been taken into account. 

The phase S(£, t) of the soliton takes the form: 

S(£, t) = -et±mu£Tmu'j- tanh" 1 [c tanh(6f)] , (5.3.61) 

with: 

c =( i+ j^r- <5 - 3 - 62) 
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a dimensionless constant. 

The EIP effect modifies the phase of the soliton. In the limit k — > 0, i.e. when 
the EIP is switched off, the phase of the soliton becomes equal to the phase of the 
soliton of the cubic Schrodinger equation. 

Finally, we remark that in the case of a pure exclusion principle (k < 0), the 
soliton exists, as we can see from (5.3.61), only if: 

Ah 2 > \n\pmN 2 . (5.3.63) 

If we take into account the maximum value of the quantity p(£), that is p(0) = 
pm N 2 /4h 2 and the maximum number of particles that can be put in a site: 

Pmax = rr , (5.3.64) 

\K\ 

Eq. (5.3.63) can be written in the form: 

p(0) < p max • (5.3.65) 

This imposes no violation of the exclusion principle in the central site, where the 
maximal occupation exists and, consequently, no violation of the exclusion principle 
on all the other points of the space. 

Taking into account Eqs. (5.2.55) and (5.3.57), we can write the potential U(p), 
which generates the soliton given by Eqs. (5.3.59) and (5.3.61), as: 

U(p) = -^p* -K-mu 2 —^- . (5.3.66) 
2 2 1 ~\- k p 
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Chapter VI 



EIP-Gauged Schrodinger Model: 
Chern-Simon Vortices 

In this chapter we describe a possible application of the model introduced in the 
chapter IV, describing systems of interacting particles obeying to EIP. Of course, 
the model can be used whenever the physical circumstances are such that collective 
excitations of charge particles occur like, for instance, in the study of degenerate 
plasmas. Here we study in some detail the properties of static, self-dual, Chern- 
Simons (CS) vortices. It was emphasized by several authors that CS theories could 
describe effects observed in the recently discovered high-T c superconductors. 
Before discussing the ElP-vortex solutions, it might be worthwhile to mention how 
such solutions were historically discovered. 

One of the most discussed topics in condensed matter is undoubtedly the supercon- 
ductivity phenomenon. In the original Ginzburg-Landau model [99], the order pa- 
rameter described by the field interacts with a Maxwell like gauge field. Although 
the equations depend on three free parameters, two of them can be eliminated by 
an appropriated scaling, leaving only one relevant physical parameter A. Supercon- 
ductors of type I and II correspond respectively to the value of the parameter A less 
or greater than one. Subsequently it was shown that the Ginzburg-Landau model 
admits vortex like solutions [125]: localized flux tube surrounded by a circulating su- 
percurrent. These vortices were experimentally observed in superconductors of type 
II. In the Bogomol'nyi limit, A = 1 vortex like solutions acquire interesting prop- 
erties, in particular static solutions are admitted because of the absence of forces 
exchanged among vortices [126]. 

Nielsen and Olesen [127] rediscovered these solutions in the context of the relativistic 
generalization of the Ginzburg-Laudau model, known as the Abelian Higgs model 
[128, 129, 130]. These authors were looking for string-like objects in relativistic field 
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theory. It turns out that these vortices have finite energy per unit length in 3+1 
dimensions (i.e. finite energy in 2+1 dimensions as the vortex dynamics is essen- 
tially confined to the x — y plane) quantized flux, but are electrically neutral and 
have zero angular momentum. In particle physics theories these solutions may be 
interpreted as strings joining confined quarks, while in cosmology theories may be 
interpreted as cosmic strings produced at a phase transition in the early history of 
our universe [131]. 

Subsequently, Julia and Zee [132] showed that the SO (3) Georgi-Glashow model, 
which admits t'Hooft-Polyakov monopole solutions, also admits its charged general- 
ization named dyon solutions with finite energy and finite, non zero, electric charge. 
It was then natural for them to enquire whether the Abelian Higgs model, which 
admits neutral vortex solutions with finite energy (in 2+1 dimensions), also admits 
its charged generalization or not. In the same paper, Julia and Zee discussed this 
question and showed that the answer is negative, i.e. unlike the monopole case, 
the Abelian Higgs model does not admit charged vortices with finite energy and 
finite and non zero electric charge. More than ten years later, Paul and Khare 
[133] showed that the Julia-Zee negative result can be overcome if one adds the CS 
term to the Abelian Higgs model. In particular, was showed that the Abelian Higgs 
model with CS term in 2+1 dimensions admits charged vortex solutions of finite 
energy and quantized, finite, Nother charge as well as flux. As an extra bonus, it 
was found that these vortices also have non zero, finite angular momentum that is 
in general fractional. This strongly suggested that these charged vortices could in 
fact be charged anyons, as it was rigorously shown by Frohlich and Machetti [134]. 
In the last years, Jackiw and Pi [115], studying the nonrelativistic reduction of the 
Abelian Higgs-CS model to the Schrodinger-CS model, where a nonlinear potential 
U(p) oc p 2 appears, were able to resolve it in the self-dual limit. In particular, 
this model admits nontopological vortex solutions in which the electric charge, the 
magnetic flux and the angular momentum are altogether quantized quantities, while 
energy and momentum are equal to zero, a feature of the self-dual solutions. The 
same model was studied in presence of an external magnetic field [135], the impor- 
tance of it being evident if we take into account the applications of the CS theory 
to the fractional quantum Hall effect [103]. 

Also the CS-Abelian Higgs model admits self-dual solutions when the nonlinear po- 
tential takes the form: U(p) oc p (p — v 2 ) 2 . This potential introduces a spontaneous 
symmetry-breaking mechanism and consequently the topological vortex makes its 
appearance in the symmetry-breaking phase. The solutions of the relativistic model 
may be considered as the analogous in (2+1) dimensions of the magnetic monopoles 
of 'tHooft-Polyakov [51, 128, 136]. 

Self-dual planar solitons can be also found in many theories of fermions where the 
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dynamics is described both in the frame of the Dirac equation with gauge field inter- 
acting by means of CS [137] or Maxwell-CS terms [138], and in the nonrelativistic 
model of the Levy-Leblond-CS equation [139]. In these theories, as in the scalar 
ones, after a calibration of the nonlinear potential, it is possible to show that each 
component of the spinorial field satisfies the Bogomol'nyi equation reducible to the 
Liouville differential equation whose solutions are known. 

An interesting generalization of the Abelian Higgs-CS model was formulated by 
Torres [140] and extended by other authors [141, 142, 143]. In this model a matter 
field is coupled in a non minimal way with a gauge field through the presence of a 
term that behaves as an anomalous magnetic moment. In the Bogomol'nyi limit the 
model admits topological and nontopological vortex solutions with non quantized 
electric charge, magnetic flux and angular momentum while in the topological sector 
the energy is quantized and infinitely degenerate. 
For an extensive review on CS theories see Ref. [144]. 



6.1 Static solutions 

Let us describe the model studied in this chapter. We consider the particular situ- 
ation in absence of the external potential V(t, x) and when the interaction of the 
gauge field is described exclusively by means of the CS term. This approximation is 
specially relevant in the context of condensed matter systems since in the long wave 
length limit, (low energy domain), the high-derivative Maxwell terms are dominated 
by the first order derivative CS one [51, 115, 136]. Thus we set 7 = and V = 
so that the motion equations for the matter field obtained from (4.1.3) and (4.1.5) 
become: 

h 2 

ihcD ^j = -—D 2 i> + A(p, J)ij + F(p)ij , (6.1.1) 

with F(p) = dU(p)/p and the nonlinear term A(p, J) given by Eq. (4.1.15). As we 
will show, the arbitrary nonlinear potential U (p) can be selected in order to permit 
the existence of self-dual vortex solutions for Eq. (6.1.1). In the self-dual limit we 
are able to decouple the gauge fields equations from the matter ones and reduce the 
evolution equation of the ijj field to an ordinary differential equation which will be 
solved numerically by means of the Runge-Kutta algorithm. 

The motion equations for the gauge fields, when only the CS term is present, become: 

g -e^F m = e -r . (6.1.2) 
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In section 4.3 we have seen that Eqs. (6.1.2) can be solved allowing us to write the 
gauge fields as a function of the sources p and J. Therefore Eq. (6.1.1) can be seen 
as an highly nonlinear Schrodinger equation for the only field ip. The model that 
we come to study is a continuous deformation, in the parameter k, of the model 
presented and studied in Ref. [51] by Jackiw and Pi. 

To search the static solutions of Eqs. (6.1.1) and (6.1.2) we make use of the 
property of invariance under gauge transformation of the Lagrangian: 



A, 



^—i(e/h c) ui 



IK 



(6.1.3) 
(6.1.4) 



We choose to work in the London gauge where the matter field ip becomes real, 
because its phase Arg(ip) is absorbed by the field A^. By using Eqs. (6.1.3) and 
(6.1.4), and setting the parameter uj = (fic/e) Arg(ip), we can see that after the 
introduction of the new fields: 



<Kx) 



IX) 



he . 



Xp(x) = A^x) + —d^Arg [if;(x)] , 



the expression of the Hamiltonian becomes: 



H 



I 



h 2 (Vp) 5 



8m 



+ 



2mc 2 



p{l + K P ) x 2 + U(p) 



d 2 x , 



(6.1.5) 
(6.1.6) 



(6.1.7) 



where p = (ft 2 is the particle density. We observe that for k < we have < p < 1/\k\ 
while for k > we have p > so that 1 + up > 0. Therefore, from (6.1.7) we see 
that, if U(p) = 0, H is a semidefinite positive quantity limited below by zero. If we 
take into account the relation: 

a±log , + 8 |£ (1 + K ^, ± |^(^) 2 + (|£) 2 (1 + Kp)x2 



2 e 

+i— (1 + Kp) 1/2 (x± d T logp - x T d ± lo gp) , 
with d± — di ± i 82 and x± — Xi =t i X2, the Hamiltonian becomes: 



(6.1.8) 



H = ^— [ <9±log p + i^ (1 + npY I 2 X ± 2 pd 2 x 

OTYl J llC 



e h 
3mc 



_ e « [ 

K J 



Xidj \(l + k p) 3/2 - l] d 2 x + 



ju(p) 



d 2 x 



(6.1.9) 
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The second term in Eq. (6.1.9) is obtained by taking into account the identity: 



(l + M V2 ^p= —d, t \(l + K pf 2 -l 

O K 1 



(6.1.10) 



This term can be integrated by part. If we neglect the surface terms and take into 
account the Gauss law given by the time component of Eq. (6.1.2): 



-gB = ep , 



(6.1.11) 



the Hamiltonian assumes the form: 

h 2 



H 



8m 
+ 



I 



2 e ,2 
d±log p + i — (I + np) 1/2 x± pd 2 x 

he 



U(p)T 



e 2 h 



3mcg k 



[l + npf 2 -l] p \ d 2 x . 



We introduce now the quantities: 



7~Co 

K(p) 



8m 
1 

K 



2e 



d± log p + i ^ ( 1 + k p) 1/2 x± 
n c 



[l + npf 2 -! 



p , 



a 



=F 



heg 



(6.1.12) 

(6.1.13) 
(6.1.14) 
(6.1.15) 



and observe that K{p) > 0, V k G M and p > 0. The Hamiltonian assumes the form: 



H = Jn d 2 x + J 



U(P) 



ti 2 a 
3m 



K{p) 



d x 



(6.1.16) 



We observe that, independently on the sign of a (which will be discuss below), the 
quantity H + (h 2 a/3 to) K(p) > 0. 

Due to the canonicity of the system, the motion equation (6.1.1) in the Hamilton 
formalism and in the case of static solutions becomes: 



5H 



which can be written as: 
5 



I 



Ho d 2 x + 



d 



Sip* J dp 



0, 



h a 



V> = o . 



(6.1.17) 



(6.1.18) 
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Following Ref. [51, 115] in which the Jackiw and Pi model (k = 0) is discussed, it 
is easy to realize that the relations: 

d± logp + (1 + k p) 1/2 X ± = ■ (6.1.19) 

U(p) + ^K(p) = 0, (6.1.20) 

ensure that the energy of the system vanishes. Eqs. (6.1.19) are two differential 
equations of the first order that permit to write the gauge field x+ an d X- as func- 
tions of the field p, while the (6.1.20) fixes the analytical expression of the nonlinear 
potential U(p). In the limit k — > 0, Eq. (6.1.19) reduces to the Bogomol'nyi equa- 
tion [126], while in the same limit the potential U(p) given by (6.1.20) becomes 
proportional to p 2 , in agreement with the results of Ref. [51, 115]. Making use of 
Eqs. (6.1.19) we can determine the field x if the expression of the field p is known: 

If we take into account the expression of the current (4.1.11), in the London gauge: 

J=-—p(l + K P )x, (6-1.22) 
mc 

and using Eq. (6.1.21), we can write J as a function of the field p: 

h 

J = T— (1 + /tp) 1/2 V Ap . (6.1.23) 
2 m 

It is easy to verify that J can be rewritten as a curl and therefore it is a fully 
transverse current. Finally, after integration of the relation: 

di A = -E i = €ij J j , (6.1.24) 



we obtain: 



h 2 a 1 



A, = Xo = o k (p)- (6- 1 - 25 ) 

6me p 

Equations (6.1.21) and (6.1.25) allow us to obtain the gauge field Xn when the shape 
of the matter field p is known. In the limit k — > 0, the field Xn assumes the expression 
given in Ref. [51, 115], so we conclude that the model described by the Lagrangians 
(4.1.3) and (4.1.5) generalizes the Jackiw and Pi model when the system obeys to 
EIP 
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Now we calculate the field p. Taking into account the relation B 
(6.1.21) it follows that: 



^ ft c„ 
B = ±— V A 

2e 



V A log p 



VA%, from 



(6.1.26) 



By using the Gauss law (6.1.11) 
equation for the field p: 



we obtain the following second order differential 



A log 



1/2 



;i + / t p) i /2 + i 



= — 2 a p , 



(6.1.27) 



which, in the limit of k — > 0, reduces to the Liouville differential equation: A log p = 
—2a p. We are not able to find the analytical solutions of Eq. (6.1.27). Numerical 
radially symmetrical solutions for the fields p will be discussed in section 6.3. Here 
we only anticipate that, how in the limit case k — 0, nonsingular, nonnegative 
solutions will be obtained when the numerical constant a is positive. Hence the =F 
sign must be chosen according that of g. 

We consider now Eq. (6.1.1) in the London gauge for the static configurations: 



_^_ApV2 
2m p 



1/2 



2mc 2 



[1 + 2k P ) X 2 + e X o + F{p) 



(6.1.28) 



It is easy to verify that the fields %, Xo an d p given by Eqs. (6.1.21), (6.1.25) and 
(6.1.27), are the required solutions. 

We remark that the self-dual character of Jackiw and Pi static solutions is not 
altered by the presence of EIP. This self-duality can be easily recognized if we define 
the following transformation: 



ip — > rj — R(r 
that changes the amplitude of the field ip in 



1/2 e -i(e/hc)S 



np) 



1/2 



4 (1 

p(r) -> R(r) = - )- 

' K (l + /tp) 1 /2 + l 

Now Eq. (6.1.19) in the 77 field becomes: 

D± V = 0. 



(6.1.29) 



(6.1.30) 



(6.1.31) 



As a consequence, the solutions of Eq. (6.1.19) are the ground states of the sys- 
tem. This can be easily demonstrated by computing the component of the energy- 
momentum tensor. We have previously noted that the energy vanishes for the 
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solutions of Eq. (6.1.19) well-behaved at infinity. The density of momentum T° l 
obtained in section 4.3, can be written by using (6.1.19) as: 



K 



1/2 



(6.1.32) 



that appears to be transverse. In the same way we can show that the flux density 
of energy and of linear momentum are: 

2 



rpiO ^ 

rpij 



e h 



2g \3mck / 
0, 



dj 



1 + npf -2(1 + npf !/2 + 1 



(6.1.33) 
(6.1.34) 



and (6.1.33) appears to be a transverse quantity. Then, apart from total derivative 
terms, the energy-momentum tensor vanishes for the solutions of Eq. (6.1.19). 
Finally, the expression of the density of angular momentum is given by: 



h c 



2hc 



e ljX *T^ = ±— V ■{[{! + K P ) 1 l 2 -l}x} T — \(l + Kp) l ' 2 -l\ , (6.1.35) 

and, if the field p vanishes at infinity more rapidly than 1/x 2 , the expression of the 
angular momentum becomes: 

2h 



< M >= T — f[(l + ^p) 1/2 -l}d 2 x 

K J 



(6.1.36) 



An alternative expression, useful in the next section, can be obtained starting from 
Eq. (4.3.40) that we rewrite for convenience: 

< M >= ~y J x A W - i> W*) d 2 x - - J x A Apd 2 x , (6.1.37) 

that, in the London gauge, assumes the form: 

< M >— — - J {x A x) pd 2 x . (6.1.38) 

Using Eq. (6.1.11) it becomes: 

< M >= ^J(x A X ) ■ (V A X )d 2 x , (6.1.39) 

which, with easy algebraic computation and taking into account the transversally 
of the field %, can be written as: 



< M >= - 

c 



2 (X 2 ) x ~ ( x ■ X) X 



A dl , 



(6.1.40) 



where the integral is to be taken both around a circle at infinity and around in- 
finitesimal contours surrounding the poles of x (zeros of p). 
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6.2 Vortex like solutions 



We study the solutions with angular symmetry (vortices) of the system. The wave 
function ip for a vortex takes the form: 



if)(r, 6) = p(r) 1 ^ 2 exp 



(6.2.1) 



where p(r) satisfies the equation (6.1.27) that in polar coordinates becomes: 

1 



d_ 

dr 



d , 4 (1 + ztp) 1 / 2 
r — log ■ 



dr 



K 



(l + /cp )l/2 + 1 



= —2a pr . 



(6.2.2) 



Analytical solutions of (6.2.2) are not known, however in the limit of k p — > Eq. 
(6.2.2) becomes the Liouville differential equation whose solutions are known [145]. 
In section 6.1 we have shown that the value of the main physical observable as- 
sociated to the solutions of Eq. (6.1.19) is zero. This statement is true with the 
exception of the mass, of the angular momentum and of the electric charge that, for 
the vortex solutions, is given by: 



Q = 2 7r e J p{r) r dr . 



(6.2.3) 



Using Eq. (6.2.2), we can write: 



Q = - 

a 



lim r f(np) — lim r finp) 



with: 



d 



log(/cp) 



(6.2.4) 



(6.2.5) 



(l + /tp) 1 /2rf r 
where only the asymptotic behavior of p is invoked. 

We remark that finite energy solutions of the system described by Hamiltonian 
(6.1.7) with the nonlinear potential given by (6.1.20) must vanish to infinity, while 
around the origin become a constant which can also be zero. We can expand the 
solution of Eq. (6.2.2) in correspondence of the zeros of p in a power series of 
p and, taking only the first order terms, we can see that the following equation 
p" — (p') 2 /p + p'/r = 0, where the prime indicates a derivative with respect to the 
r variable, must be satisfied. This equation admits power like solutions, so that we 
can write the following asymptotic behaviors of the solutions of Eq. (6.1.27): 



c or t 



p>o, 







(6.2.6) 



92 



P = C 00 r~\ 7 >0, r^oo, (6.2.7) 

where Co, Coo, P and 7 are integration constants. 

Inserting (6.2.6) and (6.2.7) into Eq. (6.2.4) we obtain the expression of the charge: 

Q = — (/? + 7 ), (6.2.8) 
a 

holding for solutions vanishing at the origin. 

Differently, we note from (6.2.4) that if the solutions p(r) at the origin are equal to 
a constant, only their behavior at the infinity determines the electric charge. In this 
case the expression of the charge Q becomes: 

Q = — 7. (6.2.9) 

a 

By comparing Eqs. (6.2.8) and (6.2.9) we can conclude that the expression of the 
charge Q is given by (6.2.8) with 7 > and f3 > 0. Taking into account Eq. (6.1.21), 
we can deduce the asymptotic behavior for r — > 00 of the gauge field x % : 

X 1 ~ ±|^7^, (6.2.10) 



r^oo 



2e ' r 



and 



X 1 ~ T^P^, (6.2.11) 
r^o 2 e r 2 

for r — > when p — > 0. The singularity of x* for r — + is due to the particular 
choice of the gauge and can be eliminated by making a gauge transformation U = 
exp(ieu/hc) with: 

u = T—PO ; 9 = arctan^ 2 /^ 1 ) • (6.2.12) 

By using the (6.2.6) and (6.2.12) we can write the expression of the field ip for small 
values of r as: 

iP(r,9) = C r /3 e ±i % e . (6.2.13) 

In order to derive the field ip as a monodrome function, the positive number /3 should 
be an even integer: 

P = 2(v-1) , (6.2.14) 

where 1/ G IV is the vorticity of the system. Therefore, we can write the charge Q 
in the form: 

Q = — [2u + 1 -2}. (6.2.15) 



a 
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In section 6.2 we have remarked that when k — > the model studied in this paper 
reduces to the model of Jackiw and Pi. Therefore, in this limit, we have 7—^2 (V+l) 
and, consequently, the discretization of the charge Q = Anev/a [51, 115]. This 
observation enables us to write 7 in the form 7 = 2($ )VK + 1) with £ )vk — v when 
k, — > 0. Therefore the electric charge Q can be written as: 

a 

In the case k ^ the parameter £ >VK is a continuous function of the boundary 
conditions. As a consequence the charge Q loses its discretization. The behavior of 
the function £ lK for a system with vorticity v — 1 will be studied in the next section. 
Finally, an alternative expression of the charge Q can be obtained from the integral 
form of Eq. (6.1.11): 

Q = -gjvA X d 2 x = -gf X -dl, (6.2.17) 

where the integral is performed both on the boundary at infinity and on the in- 
finitesimal circles around the poles of the field x (zeros of p). By using (6.2.10) and 
(6.2.11) we obtain again the (6.2.8) and (6.2.9). This result means that the poles of 
the field x (zeros of p) are placed at the origin and at the infinity. 
Finally, from Eq. (6.1.40), taking into account the asymptotic behavior of X and 
using (6.2.9) we obtain the angular momentum for the 1-vortex solutions: 

O 2 

< M >= -P— . (6.2.18) 
Aixcg 

Analogously by using (6.2.15) we obtain the following expression of < M > for the 
z/- vortex solution: 

eQ(£ +1) Q 2 
< M >= - ■ ; - — . (6.2.19) 
cga Ancg 

We observe that Eq. (6.2.19) in the case v — 1 reproduces exactly the (6.2.18) and 
then it holds for V v e IN. We recall now that the spin of the v- vortex solutions is 
given by < S >= —Q 2 /4-Kcg. From (6.2.19) we have immediately the following 
expression for the orbital angular momentum < L >=< M > — < S >: 

< L >= eg(C " + 1) . (6.2.20) 

cga 
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6.3 Numerical analysis 

In section 3.2, we have shown that the transformation: 

K — > A K , 

ip(t, x) -> Xi/;(t, x) , (6.3.1) 

can be used to rescale, in the evolution equation of the system, the value of k to 
±1. Introducing now the adimensional variables: 



2a 

n = \k\ p , (6.3.2) 

Eq. (6.2.2) becomes: 

did [l + an (y)]V2-l ] 

— \y — log 7— — / mi/2 i i r = - n (v)y » 6 - 3 - 3 

where a — k/\k\ takes the value +1 when the inclusion principle holds (k > 0, n > 0) 
and, analogously for the exclusion principle (k < 0, < n < 1), we have a = —1. If 
we take into account Eq. (6.3.2), the expression for the matter field becomes: 



fc(r) = R B (vW r ) ' (6 ' 3 ' 4) 

In order to integrate numerically Eq. (6.3.3), let us introduce the auxiliary field 

d [l + an(y)]^-l 
dy [1 + (jn(y)J 1 /^ + 1 

so that the second order differential equation (6.3.3) can be transformed into a 
system of two first order differential equations: 

^l = n(y)[l + an(y)]^ Z M^ (6 . 3 . 6 ) 

dy y 

dz(y) 



dy 



= -yn(y) ■ (6.3.7) 



The system (6.3.6)-(6.3.7) can be integrated numerically by using of the Runge- 
Kutta method, obtaining the shape of the field n(y). 
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In the following we discuss a few numerical results. In figure 6.1 are plotted the 
normalized shapes of a 1-vortex [y = 1) for the three = 0, ±1. Because of 

the proportionality between p and the magnetic field B, figure 6.1 reproduces also 
the behavior of the field B having the same polarity in each point of the space and 
a toroidal configuration around the core of the vortex. This is strictly true when 
the CS coupling constant g is negative. For positive values of g the behavior of the 
magnetic field B have opposed sign to the field p and therefore the magnetic flux 
attached to each particle is opposed to the case g < 0. This is in agreement with 
the following symmetry that hold in presence of the CS term: g — > —g, x 1 <-> x 2 , 
.1' - - A 2 . 

In figure 6.2, the electric field E r for the same 1-vortex of figure 6.1, is plotted as a 
function of r. The polarity of E r for the nontopological 1-vortex results to be always 
positive and directed radially. 

In figure 6.3 and 6.4 the shape of the fields p and E r , in the case of a 2- vortex [y = 2) 
are plotted. Now the field B has an annular shape around the core of the vortex 
and vanishes at the origin and the infinity. As a consequence, the electric field E r 
takes a polarity inversion corresponding to the points of maximum of the field B 
and it is distributed in two annular concentric regions of opposite polarity placed 
around the core of the vortex. All the curves show the effect of the EIP consisting 
in a main localization (k > 0) or derealization (k < 0) of the vortex. The same 
effect can be observed also in the shape of the electric field. Moreover, the intensity 
of the maximum of the field E r is emphasized in the attractive system respect to 
the repulsive one, in agreement to the inclusion or exclusion effect produced by the 
presence of EIP. Finally, in figure 6.5 it is shown the behavior of the function £ vk 
for an 1-vortex in the three cases k = 0, ±1 as a function of the value of the field 
p(r) at the origin, here indicated as p max , which we assume as an initial condition to 
integrate Eq. (6.3.3) (the other condition is given by dp/dr = for r = 0). We note 
that in the case k = the function £ 10 does not depend on p max and its value is 1. 
This result, obtained numerically, is equal to the one derived analytically within the 
model of Jackiw and Pi, and implies the discretization of the electric charge. On the 
contrary, in the case k — ±1 the quantity £ 10 is a continuous increasing (y = — 1) 
or decreasing (v — 1) monotonic function of p max . 

This dependence of £ x ±1 on the initial condition implies that now Q and < M > 
lose their discretization and become continuous quantities. In fact, the value of the 
charge and of the other quantities depend from the asymptotic behavior of the fields 
around the zeros of p which turn out to be function of the integration constant. 
Now, the requirement that the field ip should be a single value forces the constant 
Co to be an integer, but no condition are met for C^. 

In presence of the EIP we can observe a new fact. From figure 6.5 we note that when 
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k = — 1 the parameter and than also Q and < M > have an upper bound. So, as 
for other nonanalitically integrable models we have obtained continuous quantities 
but for a repulsive system these quantities can run in a limited range. Numerically 
we obtain £ max ~ 1.156. 

In conclusion, we have studied the static solutions of a model describing a many 
body system in the mean field approximation, obeying to a generalized exclusion- 
inclusion principle and in the presence of the Chern-Simons interaction. By selecting 
the nonlinear potential U(p) in the form (6.1.20) we have shown the existence of 
self-dual static solutions, satisfying a nonlinear first order differential equation a 
la Bogomol'nyi. The solutions are states with zero energy and linear momentum. 
Subsequently, we have considered the subset of the vortex-like solutions, obtaining 
the expression of the electric charge and the angular momentum, while their shape 
has been determinate by numerical integration of Eq. (6.3.3). The model here 
studied can be considered as a continuum deformation of the Jackiw and Pi one [51, 
115] performed by the parameter k introduced by the exclusion-inclusion principle. 
In the model of Jackiw and Pi the vortex solution takes discrete values for the 
charge, magnetic flux and angular momentum which are proportional to the vorticity 
number. In the present model these quantities take values in the continuum and 
when the system obeys to an exclusion principle (k < 0) it has an upper bound. 
Finally, we remark that in the model of Jackiw and Pi, dynamical solutions can 
be obtained from the static ones by a boosting, as a consequence of its invariance 
respect to the Galilei symmetry. On the other hand, in the frame of our model, as 
we have shown in chapter III, the EIP potential breaks this symmetry, so that we 
must specifically study nonstatic solutions. 
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Chapter VII 
Conclusions 



In this work a nonlinear canonical theory describing systems of identical particles 
obeying to a generalized exclusion-inclusion principle (EIP) has been developed. In 
the mean field approximation EIP takes into account collective effects of repulsive 
(exclusion) or attractive (inclusion) character. 

The system here considered is described by a nonlinear Schrodinger equation ob- 
tained, in the picture of the canonical quantization, from a classical model obeying 
to EIP, requiring that the quantum current density in the continuity equation takes 
the expression j K = j (1 + k p) where j is the standard quantum current density 
of the linear theory. 

To this purpose we have introduced in the evolution equation a complex nonlinear 
term A(p, j) deriving from a nonlinear potential U(if), if}*) which simulates a collec- 
tive effect among particles. The expression of this nonlinearity is strongly affected 
by the quantization method adopted. In fact, the kinetic approach here used fixes 
only the imaginary part of A(p, j). The real part is determined by the requirement 
that the system be canonical, which leads to an expression for ReA(p, j) defined 
modulo of a quantity obtainable from an arbitrary real potential U(p). This arbi- 
trarily allows to consider within EIP other interactions acting among the particles. 
After the introduction of the Lagrangian and Hamiltonian functions, we have taken 
into account the Ehrenfest relations. From their analysis we obtain that the canon- 
ical systems obeying EIP, also in the presence of an arbitrary potential U (p) and in 
the absence of external forces, are non dissipative processes (E, < P > and < M > 
are conserved) and obey to a nonlinear kinetic (d < x c > jdt — j j (1 + up) dx 
where < x c > is the mean value of the central mass of the system). Thus, EIP 
potential Ueip does not modify the dynamic behavior which can be affected only 
by the presence of an external potential V . Conversely, £7eip is responsible for the 
formation of localized stationary states (solitons). 
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These properties of the nonlinear EIP potential are obtained rigorously studying the 
symmetries of Schrodinger equation. By studying the space-time translations and 
using the Nother theorem, the expression of the energy-momentum tensor and the 
related conserved quantities have been obtained. The main results are: 

1 . EIP changes the expression of the energy density T 00 and of the fluxes T l ° and 
T l K The energy density is a semi definite positive quantity both for positive 
and negative value of k, consistently with a nonrelativistic theory. 

2. The expression of the momentum density is different from that of the current 
(P 7^ mj) and, as a consequence, the Galilei invariance is lost. 

3. The presence of U^ip introduces in the system a dimensional coupling constant 
([k] = L~ D ) and, as a consequence, the conformal symmetry of the system is 
lost. 

4. We have found a scaling transformation which permits us to reduce the cou- 
pling constant k to 1 when the inclusion principle holds, or to -1 when the 
exclusion principle holds. 

5. The discretized symmetries P and T are not lost in the presence of U E ip- 

A powerful issue, for its generality, is obtained introducing a class of nonlinear uni- 
tary transformations. In the Schrodinger equation the nonlinear term introduced 
by EIP is a complex one. By means of an appropriate transformation, it is possible 
to make real this quantity. As a consequence, the transformed system, which in 
general is described again by a nonlinear Schrodinger equation, obeys to a linear 
continuity equation where the quantum current is the same as in the linear theory. 
This method, introduced by us for EIP systems, can be applied to a large class of 
nonlinear Schrodinger equations obtained from a variational principle. 
As an application of the model describing the dynamics of a collective neutral quan- 
tum particles obeying to EIP we have studied the property of solitary wave solutions. 
These results can be used to study Bose-Einstein condensates recently obtained at 
low temperature in rarefied vapor of metals like 7 Li, 23 Na and 87 Rb. Typically these 
structures are studied by means of the cubic nonlinear Schrodinger equation. In 
place of this nonlinearity, other nonlinear terms can be considered like, for instance, 
those introduced by EIP with positive value of the k parameter to take into account 
the attractive effects due to the statistical interaction. 

We have studied systems subjected to EIP where the matter field is coupled to a 
gauge field of the abelian group U(l). 

We have also studied a theoretical model in the presence of gauge fields with a 
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dynamics described by a Maxwell-Chern-Simons Lagrangian. The most important 
properties of the system and its symmetries have been studied. We have shown that 
the anyonic behavior due to the Chern-Simons term is not disturbed by the pres- 
ence of EIP. Thus the system describes anyons obeying to EIP in the configuration 
space. 

Applications of this model can be found in the physics of condensed matter and 
in the high T c superconductors models. We have investigated stationary solutions 
when the dynamics of the gauge fields is described by the Chern-Simons term alone. 
In presence of the potential £/eip and of a nonlinear potential, real function of the 
field p, we have found self-dual solutions describing iV-vortex systems. 
The principal properties of these excited states are analytically studied. After a 
numerical analysis, we can draw the shape of the vortices and compare them to 
the ones already known in literature and deduced without considering EIP. The 
conclusive results are: 

1. The values of the charge and angular momentum of the N- vortices are con- 
tinuous, while in absence of EIP are discrete. 

2. In case of systems obeying to an exclusion principle the value of these observ- 
ables hase an upper bound and the limits can be calculated numerically. 

The numerical analysis permits us to relate the vortex profiles of systems with 
both the exclusive and inclusive effects of different weight showing how the conden- 
sate amplitude changes when the intensity of EIP changes. 
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Figure 1.1: The exchange of two particle is realized by means of rotation of one 
around the other by an angle of Atp = ± ir 



7Ti- 3 7Ti-2 TTi-i 7Ti + 7Tf + i 7Ti + +2 




7Ti-2 7Ti— l 7Ti~ 7Ti+i 7Ti+2 7T i+ 3 



Figure 1.2: One-dimensional discrete Markoffian chain where the transition proba- 
bilities are indicated. 
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Figure 5.1: Shape of soliton (in arbitrary units) with a — 0, (3 — — 1 and a = 1. 
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Figure 6.1: Plot of p for 1-vortex with k — 0, ±1 versus r [in r = (V2 a) 1 units]. 




Figure 6.2: Plot of the electric field (in Eq t = h 2 a 3 / 2 /em \pl units) for 1-vortex 
with k — 0, ±1 versus r [in ro = (V2 a) -1 units]. 
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Figure 6.3: Plot of p for 2- vortex with k — 0, ±1 versus r [in r = (V2a) 1 units]. 
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Figure 6.4: Plot of the electric field (in E 0r = h 2 a 3 / 2 /em y/2 units) for 2- vortex 
with k — 0, ±1 versus r [in r$ = (v2 a) -1 units]. 
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Figure 6.5: Plot of the parameter versus p max with k = 0, ±1. 
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